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ABSTRACT 

We surveyed the HST UV spectra of 49 low-redshift quasars for 2: < 1 C IV candidates, 
relying solely on the characteristic wavelength separation of the doublet. After consideration of 
the defining traits of C IV doublets {e.g., consistent line profiles, other associated transitions, 
etc.), we defined a sample of 38 definite (group G = 1) and five likely (G — 2) doublets with rest 
equivalent widths for both lines detected at > 'iaw^ ■ We conducted Monte-Carlo completeness 
tests to measure the unblocked redshift (Az) and co-moving pathlength (AX) over which we were 
sensitive to C IV doublets of a range of equivalent widths and column densities. The absorber line 
densityof(G = I-H2) doublets is djVciv/dX = 4.li|];^ for log A^(C+3) > 13.2, and dA^civ/dX has 



not evolved significantly since z 
of column densities /(iV(C+^)) e 



= 5. The best-fit power-law to the G 
fe(Af(C+3)/7Vo)"" has coefficient k = 



1 frequency distribution 



0.67t[J;}^ X 10- 



exponent = ~1-50^q Jgi where Nq — 10^^ cm ^. Using the power-law model of /(iV(C+'^)), 
we measured the C^'^ mass density relative to the critical density: ~ (6.20^2 52) x 10^^ for 

13 < logiV(C+^) < 15. This value is a 2.8 ±0.7 increase in r2c+3 compared to the error-weighted 
mean from several 1 < z < 5 surveys for C IV absorbers. A simple linear regression to fi^+s over 
the age of the Universe indicates that f^c+s has slowly but steadily increased from z = 5 ^ 0, 
with dflc+^/dtage = (0.42 ± 0.2) X 10"^ Gyr~\ 

Subject headings: intergalactic medium - quasars: absorption lines - techniques: spectroscopic 



1. Introduction 

With the inception of cchelle spectrometers 
on lOm-class optical telescopes, observers unex- 
pectedly discovered that a significant fraction of 
the z > 1.5 intergalactic me dium (IGM; a.k.a . 
the Lyg forest) w as enriched ( Cowie et al. 1995t 



Tvtler et al.l 119951 ). The C IV absorption lines 



have proven to be valuable transitions for study- 
ing the enrichment of the IGM since: (1) their rest 
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wavelengths AA1548. 20, 1550.77 A are redward of 
Lya A1215.67A; (2) they redshift into the optical 
for Z1548 ^ 1.5; (3) they constitute a doublet with 
characteristic rest wavelength separation (2.575 A 
or 498kms~^); and (4) have an equivalent width 
ratio 2 : 1 for Wj.,i548 ■ W^r,i550 in the unsaturated 
regime. 

Quantitative studies based on the C IV dou- 
blet measured the enrichment level to be « 10~^ 
to lO^"' the chemical abundance of the Sun 

< 



over 



the range 18 < z < 5 ([S ongaila! 12001 

Boksenberg et al.l 120031: ISchave et al.. .20031 . No 



viable model of Big Bang nucleosynthesis can ex- 
plain this observed level of enrichment; therefore, 
the metals observed have been produced in stars 
and transported to the Lya forest. The mecha- 
nisms typically invoked include "primary" enrich- 
ment by some of the earliest stars at z > 6 {e.g., 
Madau et al] I2OOII : IWise fc Abel l2008h or "con- 
temporary" injection through galactic feedback 
processes from z « 6 ^ 2 {e.g., IScannapieco et al 
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20021: lOppenheimer fc l5avil2006f) . 

Numerous high-redshift surveys have shown 
that the mass density of triply-ionized carbon rela- 
tive to the critical density of the Universe fl^+a = 
Pc+3 1 Pc.o has not evolved substantially from z = 5 
(w 1 G yr after the Big B ang) to z = 1.5 (k, 4 Gyr ; 
e.q^ Songaila _ 2001 Boksenberg et al. 2003t 



observ ations at z ~ 2.5 o ffer empirical support 



Pettini et alJ l200a iSchave et all 120031 : iSongaila 
20051) . The studies used a variety of techniques, 



from tradition a l absorption line surve ys (e.g., 
LSongaila' I2OOII: iBoksenb erg et all |2003[ ). where 
hundreds of C IV doublets were analyzed, to vari- 
ants o n the pixel optical depth (POD) method 



1 tne pixel optical deptn (rUU) metnod 
Schaveet allboOSt ISongailall2005f ). which 



statistically correlate the am ount of fl u x abs orbed 
to the C IV mass density. ISongailal (|200lf ) pio- 
neered modern fi^+s studies, and her results for 
the range 1.5 < z < 5 have been confirmed by sub- 
sequent surveys: 1.6 x 10^* < VLq+a < 3 x lO^'*. 
(We adjust all f2c+3 values quoted in this paper 
to our adopted cosmology: iJo = 70kms~^ Mpc, 
f^M — 0.3, and f^A — 0.7 and to sample doublets 
with 13 < logiV(C+3) < 15; see Appendix [C] for 
more details.) 

Recent studies have focused on increasing 
the st atistics o n C IV absorption at z > 5, 
where ISongailal (|200ll ) only detected one ab- 
sorber. These studies had to await the devel- 
opment of near-infrared spectrographs and the 
discovery of z » 6 QSOs. W i th on ly one or two 



sightli nes, iRvan- Weber et al~ ( 20061 ) and Simcoe 
( 20061 ). respectively, measured the 5.4 ^ z < 6.2 
C IV mass density to be consistent with the 
previous 1.5 < z < 5 surveys. Howev e r, the 



mos t recent work by Ryan- Weber et al.l ( 20091 ) 
and iBecker et al.l (|2009f ) observed that Oq+s at 
5-2 z < 6.2 is a factor of « 4 smaller than VIq+z 
at z < 5. Although these two studies included 
two to three times as many sightlines as Simco^ 
(|2006l) . ■ ihe results are based on < 3 detected C IV 



doublets. Obviously, the small number statistics 
at z > 5 leave the VLq+s measurements more sus- 
ceptible to cosmic variance. 

Cosmological hydrodynamic simulations have 
been used to understand the interplay between 
metallicity, feedback, the ionizin g background, 
etc., a nd t he evolution of Aguirre et al 



(|200l[ ) and ISpringel fc Hernguistl ((200^ argued 



that most of the metals observed in the IGM were 
distributed by galactic winds at 3 < z < 10, and 



Simcoe et al.l 120041) . lOppenheimer fc Dave 



(I2OO6I) evolved cosmological hydrodynamic sim 
ulations from z = 6 — > 2 with a range of pre- 
scriptions for galactic winds that enrich the IGM. 
They found that the increasing cosmic metallic- 
ity from z = 5 ^ 2 balanced the decreasing 
fraction of carbon traced by the C IV transi- 
tion (i.e., changing ionization state of the IGM). 
Thus, they neatly reproduced the observed lack of 
VLq+3 evolution. In their momentum-driven winds 
simulation (their favored 'tow' model), ^1^+3 in- 
creased from z = 6 — > 5, consisten t with the mea- 
surements available at t he time (ISongaila 2001 ; 
Rvan- Weber et all I2OO6I: ISimcod I2OO6I) and with 



the more recent re sults ( Rvan- Weber et al. 20091: 
Becker et alll2009l) . 



In lOppenheimer fc Davg (|2008l ). the authors 
included feedback from asymptotic giant branch 
(AGB) stars, in addition to a new method for de- 
riving the velocity dispersion, cr, of galaxies, which 
defined the momentum-driven wind speed. They 
compared the evolution of 17^+3 from z = 2 ^ 
in the simulations with the old and new cr-derived 
winds with AGB feedback and the new cr-derived 
winds without AGB feedback. In all three sim- 
ulations, r2c+3 did not evolve from z = 3 — > 1 
(2 Gyr to 6 Gyr after the Big Bang). In the simula- 
tion with the new cr-derived winds and AGB feed- 
back, VIq+3 increased by 70% from z = 1 ^ — 
to ~ 7 X 10~^ — over the last 8 Gyr of the 

cosmic enrichment cycle. The AGB feedback in- 
creased the star formation rate, which increased 
the mass of carbon in the IGM. This predicted in- 
crease in 51(3+3 at z < 1 can and should be tested 
by empirical observation. 

The C IV mass density at z < 1.5 has not been 
studied as extensively as at high rcdshift, where 
the C IV doublet is redshiftcd into optical pass- 
bands. At low redshift, ultraviolet spectrographs 
on space-base d telescopes are required. Recent 



studi es (e.g., iFrve et al.l [20031: iDanforth fc ShuU 



20081 ) have leveraged high-resolution, UV echelle 
spectra to examine the low-redshift IGM: the Huh- 
ble Space Telescope Space Telescope Imaging Spec- 
trograph and Goddard High-Resolution Spectro- 
graph {HST STIS and GHRS, respectively), sup- 
plemented by spectra from the Far Ultraviolet 
Spectrograph Explorer (FUSE). Through an anal- 
ysis of nine quasar sightlines observed with the 
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STIS E140M grating, iFrve et all (|2003f) found six 
C IV doublets and measured f^c+s « 12 x 10~^ 
for z < 0.1 from a preliminary analysis. With 
an expa nded survey of 28 sightlines with E140M 
spectra. lOmiforth fc ShiHll (|2008l ) detected 24 C IV 
doublets in 28 sightlines and measured flc+a = 
(7.8 ± 1.5) X 10-s for z < 0.12. These low-redshift 
flc+a measurements are consistent with that pre- 
dicted by the g-derived winds with AGB fe edback 
simulation of lOppenheimer fc Davg (|2008l ). 

These initial studies focused on C IV absorbers 
at z < 0.1 and did not present comprehensive anal- 
ysis of their survey completeness, nor did they 
take advantage of the full set of HST archival 
data. There are currently 69 sightlines in the HST 
archives with UV spectra, where the C IV doublet 
can be detected at z < 1. We have conducted a 
large survey for C IV systems in these sightlines. 
We analyzed the z < 1 data in a consistent fashion, 
which allowed for a uniform comparison through- 
out the eight-billion year interval. We introduced 
robust search algorithms and constrained the fre- 
quency distribution /(7V(C+'')) for the full sam- 
ple. Finally, we compared our results with the 
other low- and high-redshift surveys, focusing on 
the evolution of Q^+a . 

The paper is organized as follows: we present 
the spectra, the reduction procedures, and the 
measurements in § [21 our sample selection is de- 
scribed in §[3] §|l]outhnes the completeness tests; 
we analyze the frequency distribution and the C IV 
mass density in §[5] the final discussion is provided 
in § (6] and § [7] is a summary. 

2. Data, Reduction, and Measurements 

To assemble our target list, we searched the 
HST STIS and GHRS spectroscopic archives for 
objects with target descriptions including the 
terms, e.g., QSO, quasar, Seyfert, etc. Our fi- 
nal list included 69 objects with redshifts 0.001 < 
Zem < 2.8. We retrieved all available spectra 
from the Multimission Archive at Space Tele- 
scope0 (MAST), including supplementary FUSE 
data, when available (see Table [J). We reduced 



(in the case of the FUSE spectra), co-added, and 
normalized the spect ra with similar algorithms as 



described in Cooksc v et al. ( 20081 ). The reduced 



HST spectra were retrieved from MAST directly. 

One goal of this study was to search the entire 
HST archive for C IV absorption; thus, we ini- 
tially included the full archival dataset. However, 
we excluded targets that only had un-normalized 
spectra with signal-to- noise ratio S/N < 2pix^^. 
The signal-to-noise ratio was measured by fit- 
ting a Gaussian to the histogram of flux fx and 
noise cr/^ per pixel, clipping the highest and low- 
est outliers and demanding fx > 0. Ten sight- 
lines do not meet the S/N criterion: QJ0640- 
5031; TON34; HE1104-1805A; HE1122-1648; 
NGC4395; Q1331+170; IR2121-1757; HDFS- 
223338-603329-QSO; AKN564; and PG23024-029. 
They are indicated by 'Exclude (low S/N)' in Ta- 
ble[TJ Q1331-I-170 has a damped Lya sy stem, with 
molecular hydrogen lines, at z — 1.7765 ( Cui et al] 
,2005) that made the E230M spectra exceptionally 
difficult to analyzeH Therefore, Q1331-t-170 was 
also excluded, though it had S/N 5pix^^. 

In addition, we excluded the spectra for which 
higher-resolution spectra covered the same wave- 
length range (noted as 'Exclude (overlap)' in 
Table [T]). Typically, we excluded spectra with 
resolution R < 20,000 (FWHM > 15kms"^). 
There were ten targets that had no C IV cov- 
erage, given the spectral wavelength range and 
Zem- Q0026+1259, TONS180, PKS0558-504, 
PG1004-I-130, HE1029-f-140, Q1230+0947, 
PG1307-I-085, MARK290, Q1553-I-113, and 
PKS2005-489. 

Ultimately, 49 targets had spectra with usable 
wavelength coverage and S/N ratios. All reduced, 
co-added, and normalized spectra are available on- 
lineH even those not explicitly searched in this pa- 
per. The number of spectra covering the C IV 
doublet redshift range is shown schematically in 
Figure [T] 



^See htt p : / / archive . s tsci .edu/j 

^We have adopted the target name that MAST used. 
So note: B0312-770 is also PKS0312-770; QSO- 
123050+011522 is also Q1230-0115; and PHL1811 is also 
known as FJ2155-092. 



^ Q1331-|-170 also ha d a known Mg II system at z = 0.7450 
l lEllison et al.ll2003l) . which likely had blended C IV ab- 
sorption. 

■^See http://www.ucolick.org/~xavier/HSTCIV/ for the 
normalized spectra, the continuum fits, the C IV candi- 
date lists, and the Monte-Carlo completeness limits for all 
sightlines as well as the completeness test results for the 
full data sample. 
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Table 1 
OBSERVATION SUMMARY 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


Target 


RA 


Dec 




Instr. 


Grating 


R (FWHM) 


S/N 


Amin 


A„ax 




PID 


Notes 




(J2000) 


(J2000) 








(kms^^) 


(pix-i) 


(A) 


(A) 


(ks) 






MRK335 


00:06:19 


+20:12:10 


0.0258 


STIS 


E140M 


45800 (7) 


7 


1141 


1710 


17 


9802 












GHRS 


G160M 


20000 (15) 


34 


1222 


1258 


15 


3584 












FUSE 


LWRS 


20000 (15) 


< 10 


904 


1188 


99 


PlOl 




Q0026+1259 


00:29:13 


+ 13:16:04 


0.1450 


GHRS 


G270M 


20000 (15) 


5 


2785 


2831 


5 


3755 


Exelude (coverage) 










FUSE 


LWRS 


20000 (15) 


< 2 


904 


1188 


20 


Q206 




TONS180 


00:57:20 


-22:22:56 


0.0620 


STIS 


G140M 


12700 (24) 


25 


1244 


1298 


3 


7345 


Exclude (coverage) 












G230MB 


9450 (32) 


9 


2758 


2912 


1 


9128 


Exclude (coverage) 










FUSE 


LWRS 


20000 (15) 


< 5 


904 


1188 


25 


D028; PlOl 




PG0117+213 


01:20:17 


+21:33:46 


1.4930 


STIS 


E230M 


30000 (10) 


8 


2278 


3072 


42 


8673 




TONS210 


01:21:51 


-28:20:57 


0.1160 


STIS 


E140M 


45800 (7) 


6 


1141 


1710 


22 


9415 














E230M 


30000 (10) 


4 


1988 


2782 


5 


9415 














G230MB 


9450 (32) 





2759 


2913 


2 


9128 


Exclude (fi) 










FUSE 


LWRS 


20000 (15) 


< 8 


904 


1188 


53 


P107 





Note. — The targets and spectra included in the current survey. The targets, their coordinates, and their redshifts are Usted in Columns 1-4. The instrument and 
grating of the archived observations (Columns 5 and 6) are listed with the resolution R (and FWHM in kms~^), exposure time t^^p, and the proposal identification 
number (PID) from the MAST query (Columns 7, 11, and 12). The signal-to-noise ratio S/N was measured from the normalized spectra (Column 8). Columns 9 and 
10 are the wavelength coverage of each spectrum. The status of the spectra (i.e., whether it was excluded and why) is listed in Column 13. Tablef^is published in its 
entirety in the electronic edition of the Astrophysical Journal. A portion is shown here for guidance regarding its form and content. 




the-Fly Repr ocessing; and co-adde d with XIDlH 



COADSTIS (ICooksev et al.l l20nah . The STIS 



Fig. 1. — Schematic of redshift coverage for the 
current survey. The number of spectra with cov- 
erage of the C IV doublet is shown as a func- 
tion of 21548, the redshift of the 1548 hne (black 
histogram). The STIS E140M spectra covered 
^1548 < 0.1. The STIS E230M spectra typically 
covered 0.4 < Z1548 < 1. The redshift range 
0.1 < 2^1548 ^ 0.4 was covered by some E230M 
spectra as well as the medium-resolution gratings 
and GHRS spectra. The redshifts of the doublets 
detected with Wr > 3 aw^. hi both lines are shown 
with the hashes across the top. The top and mid- 
dle rows indicate the redshifts of the 27 unsat- 
urated (top) and 11 saturated (middle) doublets 
in the definite group (G = 1; see ij 13. 3p . The bot- 
tom row shows the redshift of the five unsaturated 
likely (G = 2) doublets. 

2.1. HST STIS and GHRS 

The HST spectra drive the target selection be- 
cause its UV instruments, STIS and GHRS, have 
the wavelength coverage to detect C IV in the 
z < 1 Universe. We preferred higher-resolution 
data over lower-resolution data to resolve the C IV 
doublets and better distinguish them from coinci- 
dent Lya features. The STIS echelle gratings were 
the preferred set-up in most cases, since they were 
high resolution and covered a large wavelength 
range. 

All STIS observations were reduced with Cal- 
STIS 2.23 on UT 06-October-2006 with On- 



long-slit spectra were co-added by re-binning the 
individual exposures to the same wavelength array 
and then combining the error-weighted flux. For 
the echelle data, the observations for each order 
were co-added separately, in the same manner as 
the long-slit spectra. Then all of the orders were 
co-added into a single spectrum. 

The STIS echelle gratings E140M and E230M 
contributed the most to the total pathlength of 
this survey. All 49 targets have at least one 
echelle spectrum. The E140M grating covered 
the wavelength range 1140A < A < I7IOA or 
Z1548 ^ 0.1 and had a resolution of i? = 45,000 
(FWHM = Tkms"^). The E230M grating cov- 
ered « 800 A per tilt over the range 1570 A < 
A < 3110 A or 21548 < 1. Typically, the obser- 
vations covered 2280 A < A < 3070 A or 0.5 < 
^1548 ^ 1- The E230M grating had resolution 
R = 30,000 (lOkms"^). The only other STIS 
gratings included in the C IV search were G140M 
and G230M, which have wavelength coverage simi- 
lar to E140M and E230M, respectively, albeit with 
less total coverage. For m ore information about 
STIS, see lMobasheJ ()2002f ). 

The calibrated GHRS spectra were retrieved 
from MAST. They were co-added hkc the STIS 
long-slit spectra. Observations with central wave- 
lengths that varied by less than 5% of the total 
wavelength coverage of that grating were co-added 
into a single spectrum. This avoided the problem 
of automatic scaling to the highest-S/N spectrum. 

The GHRS ECH-A, ECH-B, G160M, G200M, 
and G270M gratings were used in the final C IV 
search. The ECH-A and ECH-B gratings had res- 
olution R ^ 100,000 (FWHM = 3kms"^). The 
G160M, G200M, and G270M gratings had reso- 
lution R = 20,000 (15kms"^). We excluded all 
GHRS spectra where there existed STIS coverage 
since, in general, the STIS data have higher resolu- 
tion and larger wavelength c overage. For more in - 
formation about GHRS, see lBrandt et al.l (|l994r ). 



2.2. Far Ultraviolet Spectroscopic Explorer 

The raw FUSE observations were retrieved 
from MAST. We reduced the spectra with a mod- 



°See |http:// www.ucolick.org/~xavier/IDL/ 1 
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ified CalFUSEH v3.2 pipeline (see ICooksev et al 
20081 ). To summarize, the partially processed data 
from the exposures (i.e., intermediate data file or 
IDF) were combined before the extraction win- 
dow centroid was determined. Then, the bad pixel 
masks were generated, and the spectra were opti- 
mally extracted. 

The observations were co-added with 
FUSE_REGISTER. The FUSE wavelength solu- 
tions were not shifted to match the STIS wave- 
length solution. The FUSE observations were 
only searched for absorption lines {e.g., Ly/3, C III, 
O VI) to supplement candidate C IV systems. 
The blind doublet search made allowance for off- 
sets in redshift, due to multi-phase absorbers or 
misaligned wavelength solutions, when it assigned 
additional absorption lines to the candidate C IV 
systems (see 13. II for more details). 

We excluded FUSE spectra with S/N ratios 
too low for continuum fitting: TON28, NGC4395, 
QSO~-123050-f-011522, Q1230-t-0947, NGC5548, 
PG1444-I-407, and PG1718-h481. 

2.3. Continuum Fitting 

The continuum for each spectrum was fit 
semi-autoniatically with a B-spline. "Semi- 
automatically" refers to the subjective nature of 
continuum fitting. The best B-spline was based 
on the input parameters, primarily, the low- and 
high-sigma clips and the breakpoint spacing. How- 
ever, the authors varied these parameters until 
they agreed the "best" B-spline was a good match 
to un-absorbed regions. 

The breakpoint spacing was typically « 4 A to 
6 A. The breakpoint spacing was refined to in- 
crease in regions where the flux was not varying 
substantially and to decrease in regions where it 
was. Regions of great variation were determined 
by comparing the median flux of each bin to the 
mean variance-weighted flux of all bins, where the 
bins are determined by the initial breakpoint spac- 
ing. If the difference between the median flux in 
a bin and the mean flux were more than one stan- 
dard deviation of the flux bins, the breakpoint 
spacing was decreased. For most spectra, the low- 
and high-sigma clips were 2a and 2.5a, respec- 
tively. The regions masked out by the sigma clip- 
ping were increased by two pixels on both sides. 



The B-spline was iteratively fit to the variance- 
weighted flux in the bins, with the clipped regions 
masked out, until the percent difference across the 
continuum fit was less than 0.001% compared to 
the previous (converged). For several spectra {e.g., 
the NGC galaxies) , the semi-automatic continuum 
fit was adjusted by hand. 

As mentioned previously, the "best-fit" contin- 
uum was a subjective judgment. Several, slightly 
different continuum fits would have satisfied the 
authors. To gauge the systematic error intro- 
duced by the subjective nature of continuum fit- 
ting, we measured the differences due to chang- 
ing the continuum from the semi-automatic fit to 
one generated "by hand" for one sightline (see 



Cooksev et al. 2008f) . The root-mean-square frac- 



^ See [ftp: / / fuse.pha.jhu.edu/fuseftp/calfuse71 



tional difference in the observed equivalent width 
and column density were < 10%. 

2.4. Redshift, Equivalent Width, and Col- 
umn Densities 

We measured redshifts from the optical depth- 
weighted central wavelengths; equivalent widths 
from simple boxcar summation; and column den- 
sities from the apparent optical depth method 
(AODM; ISavaee fc Sembachlll99lh . To minimize 
the effect of spurious, outlying pixels on the mea- 
surements, the flux was trimmed: — c/^ < /a < 
1 + af^. Outlying flux pixels were set to the ap- 
propriate extrema. The trimming affected mea- 
surements in the lowest S/N regions of the spec- 
tra, typically the edges. It also assisted the iden- 
tification of absorption lines in the noisy regions 
by the automated feature search (see § 13. 1|) . from 
which the candidate C IV doublets were drawn. 
When the final C IV sample was defined and the 
wavelength bounds visually confirmed or changed, 
this flux trimming had negligible effect. For opti- 
cal depth measurements, the minimum flux was 
set, so that: fx > (0.2cr/^ > 0.05). This pre- 
vented the optical depth, r = ln(l//A), from be- 
ing overwhelmingly large for the saturated pixels. 
These cases were then reported as lower limits to 
7V(C+3). 

We measured redshifts from the mean optical 
depth- weighted central wavelengths of the absorp- 
tion lines and the rest wavelength of the transi- 
tion. The wavelengths per pixel Xi were weighted 
by their optical depth per pixel = hi{l/f\i) 
within the bounds of the absorption line, deflned 
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by the wavelength range: A; < A.; < A/i. Thus, the 
pixels with the strongest absorption dominated 
the redshift estimate: 



1 + Zab 



(1) 



The rest equivalent width Wr were measured 
with a boxcar summation over the wavelength 
bounds of the feature: 

Wr = Y.(^-fx^)6X. (2) 

V-L ~r Zabs ) , 



2 ^/i 



where 6Xi (A) is the wavelength pixel scale of the 
spectrum. The second equation is the variance of 
the Wr measurement from error propagation. The 
observed equivalent width Wobs and error crw^bs 
were measured with the previous equations but 
without the 1 + Zabs factor. 

Most of the column densities were measured 
with the apparent optica l depth method (AODM; 
Savage fc Sembachlll9911 ): 



2 

^AfAOD 



10 



14.5762 ^fe 



foscXr 



1 



J2^n( — )6v,,Ucm-') (3) 



-|^q14.5762 \ 2 A,^ 



fosc 



Sv 



SX, 



Xr{l + Zabs) J' 



where fosc (unitless) is the oscillator strength of 
the transition with rest wavelength Xr (A) and 
Sv\^i (kms^^) is the velocity pixel scale of the 
spectrum The atomic da t a and sources are tab- 
ulated in iProchaska et aL The AOD pro- 
files were used as a diagnostic (see § 13. 2p and were 
constructed from the un-summed versions of the 
above equations, smoothed over three pixels (see 
Figure [11). 

In some cases, only a column density limit could 
be set, when the AODM resulted in a measure- 
ment A'aod < '^o'Naov (resulting in an upper 



limit) or the line was saturated (lower limit). In 
low-S /N spectra, some lines satisfied both criteria 
{e.g., Z1548 = 0.40 2 2 7 doublet in the PKS0454-22 
sightline). In which case, we counted the line as an 
estimate of the upper limit if Wr < 3 aw^ and as 
a lower limit otherwise. Since the column density 
from the AODM was a poor measurement in the 
aforementioned instances, we estimated the col- 
umn density by assuming the Wr reflects the col- 
umn density from the linear portion of the curve 
of growth (COG). Then, we use whichever column 
density measurement resulted in the more extreme 
limit. For example, if the COG column density 
was lower than the AODM column density for a 
saturated, > 3aWr feature, the COG column den- 
sity was used. 

For analyses pertaining to the equivalent width, 
we used Wr, 154s from the C IV 1548 hue. For 
analyses relating to the column density A^(C+^) 
we either used the error-weighted average of A^i548 
and A^i55o when both were measurements or con- 
strained the value based on the limits of the two 
doublet lines. The greater lower limit (or the 
smaller upper limit) was used. In a few cases, 
A^i548 and iVi55o constrained a N{C'^^) range (see 
the bracketed values in Table O. For the column 
density analyses, we took the average of the limits 
and used the difference between the average and 
the values as the errors. For the 2:1548 = 0.38152 
C IV doublet in the PKS0454-22 sighthne, the two 
line limits did not overlap; we increased/decreased 
the limits by la to constrain A^(C"'"'^). 

For several sightlines {e.g., the FUSE spectra 
or the two E230M grating tilts for PG1634-f 706), 
there were overlapping spectra. We quote mea- 
surements from the spectrum where Wr was 
measured with the higher estimated significance 

Wr/(7Wr- 

3. Sample Selection 

We conducted a "blind" survey for C IV dou- 
blets, where candidate C IV absorbers and any 
associated lines were identified exclusively by 
the characteristic wavelength separation of the 
C IV doublet and the measured redshift of the 
C IV 1548 line. This eliminated any bias asso- 
ciated with identifying Lya absorbers first and 
then searching for associated C IV doublets. 
Also, for several sightlines, we did not have the 
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wavelength coverage to search for Lya. Though 
C IV sy stems frequen tly sho w strong Ly a ab- 
sorption lEUison et al...l999. : Simcoe et al.i [2004). 
C IV systems with w eak Lya absorption do occur 
(jSchave et alJuOOTD . and they might occur with 
higher frequency at z < 1 if the IGM is more 
highly enriched. Once we searched for the can- 
didate C IV doublets, we used other diagnostics 
(see § 13. 2|) and visual inspection to define our final 
sample. 

3.1. Automatic Line Detection and Blind 
Doublet Search 

First, we searched for absorption feature^ with 
observed equivalent widths W^hs > -icrw ,. wit h 
an automated procedure ( Cooksev et al] l2008[ ). 
We convolved the flux with a Gaussian with width 
equal to the FWHM of the instrument to yield 
/g- Then, adjacent convolved pixels with signif- 
icance > l<T/g were grouped into absorption fea- 
tures. The observed equivalent width for all fea- 
tures were measured with the boxcar extraction 
window defined by the wavelength limits. A; and 
Xh- The final product was a list of candidate ab- 
sorption features with Wobs > 3 <Jw„ts ' defined by 
their fiux decrement-weighted wavelength ccntroid 
and wavelength limits. 

Weighting with the flux decrement 1 — fx, where 
/a < 1, was better behaved than weighting with 
the optical depth r = ln{l/f\). For strong ab- 
sorption lines, the two weighting methods result 
in comparable values. However, for weak absorp- 
tion lines and noisy absorption features, the flux 
decrement- weighted centroid was closer to the cen- 
tral wavelength, as defined by the bounds A; < 
A < A;, 

There was no attempt to separate blends at this 
level. This method automatically recovered nearly 
all of the features we would have identified visu- 
ally, but it was not sensitive to broad , shallow fea- 
tures, i.e., b > 50kms^^. However, Ranch et al 



l|1996l ) measured the distribution of Doppler pa- 



rameters of C IV doublets at high redshift to be 
5kms~^ <b< 20kms"^. This gives us reason to 



''^We refer to absorption "features" to indicate the results 
from the automated feature-finding program, which also re- 
turned absorption line-like artifacts in the spectrum {e.g., 
regions with poor continuum fit, multiple pixels with spu- 
riously low fiux) in addition to real absorption lines. 



expect that our algorithms would miss very few (if 
any) systems. 

The automatically-detected features were paired 
into candidate C IV doublets based purely on the 
characteristic wavelength separation (2.575 A or 
498kms~^ in the rest frame). Every feature with 
Xcent between A1548 and Ai548(l -t- z^) was as- 
sumed to be C IV 1548 absorption|f| If there 
were an automatically-detected feature between 
the wavelength limits (A/,A/i) at the location of 
the C IV 1550 line, it was assumed to be C IV 
1550 absorption (see Figure [2|). Otherwise, the 
region that would have included C IV 1550 was 
used to give an estimate of the upper limit of the 
column density and equivalent width. The C IV 
1550 region was set by the wavelength bounds of 
C IV 1548: 



A/, 1550 
Aft,, 1550 



-^;,1548 
Aft ,1548 



Ar,1550 
Ar,1548 , 

Ar,1550 
A 



(4) 



r,1548 



Similarly, the wavelength bounds of the candidate 
doublet were adjusted so that they were aligned 

{e.g., W;,i548 = 'y/,155o)- 

The doublet search was then performed in re- 
verse, where automatically-detected features were 
assumed to be C IV 1550 lines, if not already in- 
cluded. The corresponding region for the C IV 
1548 absorption yielded an estimate for the upper 
limit of the column density and equivalent width. 

Automatically-detected features at the ex- 
pected locations, defined by the redshift of the 
C IV 1548 lines, of several common transitions 
were added to the candidate C IV systems. The 
common transitions were: Lya; Ly/?; C III 977; 
O VI 1031, 1037; Si II 1260; Si III 1206; Si IV 1393, 
1402; N III 989; and N V 1238, 1242. If there were 
no automatically-detected candidate Lya absorp- 
tion, the region where the Lya line would have 
been was included as an estimate of the upper 
limit of the column density and equivalent width, 
if the wavelength coverage existed. 

There were some automatically-detected fea- 
tures with Ai548 < A < Ai548(l + Zem) that were 



'Actually, we searched for doublets at higher redshifts than 
z^rn as a secondary check to the adopted values of z^rn- We 
did not find any C IV doublets with Z1548 > Zem- 



8 



-200 -100 100 200 

Relative Velocity (km s"') 
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Fig. 2. — Example of an initial candidate C IV doublet and the final Z1548 — 0.91440 system in the 
PG1630+377 sightline. Two automatically-detected absorption features were paired as a candidate C IV 
doublet based on their wavelength separation (left panel). The wavelength bounds of the candidate 1550 
line were expanded to match the range defined by the candidate 1548 line. Several automatically-detected 
features were added to the candidate C IV system as common transitions, based on their observed wave- 
lengths and the redshift of the candidate 1548 line. Finally, we visually inspected the candidate, deemed it 
a definite system, and modified the transitions and wavelength bounds (right panel). The bottom panel on 
the left shows the AOD profile of the 1548 and 1550 fines (black solid and red dotted lines, respectively; see 
13. 2p . For a description of the velocity plots, see Figure [3] 



more than 498kms~^ wide, which is the charac- 
teristic separation of the C IV doublet. For ex- 
ample, there is a damped Lya system in the spec- 
trum of PG1206+459 at Zabs = 0.92677, which 
has a strong, multi-component C IV doublet that 
is 633 km s^^ wide, and the doublet lines 1548 and 
1550 are blended with each other (see Appendix 
|A|) . These cases were visually evaluated and in- 
cluded as candidate C IV doublets, since no au- 
tomated attempt was made to separate blended 
absorption features. The search for common tran- 
sitions then proceeded as described previously. 

Finally, the redshifts, rest equivalent widths, 
column densities, and errors were measured for all 
transitions in the candidate C IV systems (see § 
12. 4[) . The measured quantities for the candidate 
doublets are listed in Tabled 



3.2. Machine-Generated Diagnostics 

We assembled several thousand candidate C IV 
systems by the process described in the previous 
section. We then used machine-generated diag- 
nostics to assist in identifying the true C IV sys- 
tems visually. The candidate systems were as- 
signed flags based on nine criteria that leveraged 
our knowledge of the C IV doublet and other de- 
sirable characteristics (see Table IS]). The crite- 
ria were the following: C IV 1548 and 1550 had 
Wr > "iow-r (flags = 256 and 128, respectively); 
the Wr ratio of 1548 to 1550 was between 1 : 1 and 
2:1, accounting for the error (64, see Equation 
[S]); the redshifts of the doublet lines were within 
lOkms"^ (32, see Equation [5]) ; there was a can- 
didate Lya absorber detected with Wr > Saw^ 
(16); C IV 1548 was outside the Lya and the H2 
forest (8 and 4, respectively); the apparent optical 
depth profiles of the doublet were in agreement 
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Table 2 

C IV CANDIDATES SUMMARY 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) (10) 


^1548 




A,. 


A, 




Wr 




\ogN 


o"log JV Flag 




(kms^^) 


(A) 


(A) 


(A) 


(mA) 


(mA) 














MRK335 (Zem 


= 0.026) 








-0.00162 




1548.19 


1545.33 


1545.97 


< 48 




< 13.09 


157 




-11.7 


1550.77 


1547.90 


1548.54 


286 


21 


> 14.39 




-0.00092 




1548.19 


1546.49 


1547.03 


153 


20 


13.72 


0.07 335 




-18.0 


1550.77 


1549.06 


1549.61 


< 59 




< 13.50 




0.00000 




1548.19 


1547.90 


1548.54 


286 


21 


> 14.09 


495 




4.2 


1550.77 


1550.47 


1551.12 


183 


31 


14.15 


0.09 








PG0117+213 (ze 


„ = 1.493) 








0.47166 




1548.19 


2278.01 


2278.98 


691 


78 


> 14.23'' 


288 




3.7 


1550.77 


2281.80 


2282.77 


< 119 




< 13.94 




0.47372 




1548.19 


2280.33 


2282.56 


< 191 




< 0.30 


128 




-23.5 


1550.77 


2284.12 


2286.36 


544 


73 


> 14.66 




0.47605 




1548.19 


2283.79 


2286.36 


522 


79 


> 14.35 


386 




-50.0 


1550.77 


2287.59 


2290.16 


978 


71 


> 15.05 










TONS210 (z™ 


= 0.116) 








-0.00073 




1548.19 


1546.84 


1547.40 


155 


20 


> 13.85 


367 




2.1 


1550.77 


1549.41 


1549.98 


< 64 




< 13.63 




-0.00001 




1548.19 


1547.64 


1548.67 


480 


25 


> 14.32 


495 




-7.9 


1550.77 


1550.21 


1551.24 


328 


36 


14.42 


0.06 


0.00090 




1548.19 


1549.43 


1549.73 


74 


23 


< 13.23 


269 




-18.0 


1550.77 


1552.01 


1552.31 


< 38 




< 0.30 





^log N measured by assuming Wr results from the linear portion of the COG. 



Note. — Summary of CIV doublet candidates by target and rcdshift of CIV 1548. Upper limits are 2(7 
limits for both Wr and log, N . The binary flag is described in Table [3] and ^ 13.21 Table [2] is published in 
its entirety in the electronic edition of the Astrophysical Journal. A portion is shown here for guidance 
regarding its form and content. 
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PG0117+213 (01:20;17.25) LyaForest Zi.^s ^ 0.57632 




-100 100 

Relative Velocity (km s"') 



Fig. 3. — Velocity plot of G = 1 C IV system in 
the PG0117-I-213 sightlinc. The regions of spec- 
tra around each absorption line are aligned in ve- 
locity space with respect to the rest wavelength 
of the transition and 2:1543 = 0.49907. Saturated 
transitions are indicated with the (red) 'S.' The 
regions used to measure Wr and logiV(C"'"^) are 
shown by the dark outline. The flux at zero and 
unity are shown with the dash-dot lines (blue and 
green, respectively); the vertical dashed line indi- 
cates V = Okms~^ corresponding to the optical 
depth- weighted velocity centroid of the C IV 1548 
transition. (The velocity plots for all G = 1 and 
G = 2 absorbers are in Appendix [Xl) 



(2); and other candidate transitions beside Lya 
were detected with W,. > 3(JWr (!)■ 

The equivalent width ratio of C IV 1548 to 1550 
was measured as follows: 



Rw — 



2 



Wr 



1548 



r,1550 



(5) 



i?2 



w 



f / 0'ty^,1548 \^ / 
W^r,1548 J V 



r,1550 



where the variance a\ 



resulted from the prop- 
agation of errors. The actual diagnostic flag = 
64 applied to doublets where 1 — or^ < Rw < 
2 -|- CTij^, . The upper limit was set by the charac- 
teristic ratio for an unsaturated doublet, and the 
expected ratio decreases for increasingly saturated 



doublets until Wr s548 — WrS550, hence, the lower 
limit. 

The redshift of the 1548 line Z1548 must equal 
the redshift of the 1550 line zisso for an un-blended 
C IV doublet, detected with infinite S/N and infi- 
nite resolution. To accommodate blending and the 
range of S/N ratio and resolution of the spectra, 
the binary diagnostic flag = 32 applied to doublets 
with |(5wciv| < lOkms^^, where: 



6vc IV 



f ^1550 ~ ^1548 
V 1 + 21548 



(6) 



A candidate system was outside the Lya for- 
est when Ai548(l + 2:1548) > Aa(l + Zem) (flag = 
8). Similarly, a candidate was outside the H2 for- 
est when Ai548(l + 21543) > Ajj^ (1 -I- Zem) (flag = 
4), where Au = 1138.867A for H2 B0-0P(7), the 



last H 2 transition considered (see iProchaska et al 



2004 for atomic and molecular line data). 

The profiles of the 1548 and 1550 lines were con- 
sidered to have matched when they agreed within 
la for 68.3% of the range A; to A,, (flag = 2). The 
profiles were defined by the apparent optical depth 
(AOD) column densities, similar to Equation [3l 



N,. 



2 



ir)14.5762 / -i 

— ^In — )5z;A,,(cm-2) (7) 

Josc^r \ JXi 



/lO 



14.5762 



\ fosc^r fxi 



■6v\ 



(cm ^) 



where i indicates the pixels between A; and Xh- 
The profiles were smoothed over three pixels to 
minimize the effects of spurious pixels and noisy 
spectra. C IV doublets with the diagnostic flag = 
2 have: 



1548 



1.1550 



< 



(8) 



for more than 68.3% of the range A/ to A^. The 
AOD profiles were computed only for the candi- 
date C IV doublets. 

The flags are binary and sum uniquely. They 
are ordered so that doublets satisfying more of the 
diagnostic criteria have higher flag values. For the 
current study, we focus on the doublets with both 
lines detected at > 3aw^, or flag > 384. 
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Table 3 
C IV DIAGNOSTIC FLAGS 



Flag 


Description 


256 


^i',1548 > 3(7Wi,,1548 


128 


Wr,1550 ^ 3crvKr,1548 


64 


1 — '^Rw — < 2 + CTRj^ 


32 


liS-ucivi < lOkms-1 


16 


Candidate Lya with 1215 > 3c7"w^,i2i5 


8 


C IV 1548 outside Lya forest 


4 


C IV 1548 outside H2 forest 


2 


> 68.3% AOD profile per element (pixel) in l-o" 




agreement 


1 


Other candidate absorption lines with Wr > 




3^ 



Note. — Binary diagnostic flags for C IV doublets. (For 
more information, sec ^ 13.21 ) 
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Table 4 
C IV SYSTEMS SUMMARY 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) (11) 


^1548 


&Vaba 












log AT 


flog N 


G Flag 




(kms-1) 


(A) 


(A) 


(A) 


(mA) 


(mA) 














PG0117+213 ( 


^em = 


1.493) 








U. 51904 




1548. 19 


2352.30 


2353.07 


94 


24 


13.47 


0.11 


2 288 




-3.5 


1550.77 


2356.21 


2356.98 


< 45 




< 13.39 






U.OYDOZ 


' ' ' 


1548.19 


2439.35 


2441.47 


728 


28 


> 14.56 


' ' ' 


1 450 




12.6 


1550.77 


2Ai3.A\ 


2445.24 


442 


25 


14.52 


0.03 










PKS0232-04 ( 


Zem — 


1.440) 








0.73910 




1548.19 


2691.50 


2693.40 


372 


34 


14.16 


0.08 


1 450 




18.7 


1550.77 


2695.98 


2697.88 


359 


41 


> 14.43 






0.86818 




1548.19 


2891.94 


2892.69 


69 


22 


< 13.13 




2 358 




-5.6 


1550.77 


2896.75 


2897.50 


< 37 




< 13.32 












PKS0405-123 (zem = 


0.573) 








0.36071 




1548.19 


2106.28 


2107.05 


86 


17 


13.37 


0.09 


1 383 




6.4 


1550.77 


2109.79 


2110.55 


< 35 




< 13.27 








53.0 


1215.67 


1653.80 


1655.51 


769 


35 


> 14.47 








15.7 


1025.72 


1395.47 


1396.23 


270 


10 


> 14.93 








8.7 


977.02 


1329.24 


1329.69 


142 


9 


> 13.68 








19.9 


1206.50 


1641.62 


1641.94 


73 


12 


12.70 


0.09 










PKS0454-22 ( 


Zem = 


534) 








0.20645 




1548.19 


1867.44 


1868.12 


142 


24 


13.64 


0.08 


2 366 




0.9 


1550.77 


1870.54 


1871.23 


81 


26 


13.70 


0.14 




0.24010 




1548.19 


1919.35 


1920.52 


644 


57 


> 14.20* 




1 494 




7.7 


1550.77 


1922.54 


1923.71 


495 


54 


> 14.39'' 






0.27797 




1548.19 


1978.18 


1979.04 


274 


52 


> 13. 83* 




1 366 




-9.6 


1550.77 


1981.47 


1982.33 


< 104 




< 14.20 






0.38152 




1548.19 


2138.54 


2139.20 


96 


32 


< 13.55 




1 494 




1.2 


1550.77 


2142.10 


2142.76 


110 


27 


> 13. 73* 








3.8 


1215.67 


1679.22 


1679.74 


< 466 




< 14.23'* 






0.40227 




1548.19 


2170.55 


2171.37 


177 


29 


< 13.72 




1 494 




1.7 


1550.77 


2174.16 


2174.99 


140 


31 


> 13.84'' 








-5.4 


1215.67 


1704.35 


1704.99 


< 458 




< 14.23* 






0.42955 




1548.19 


2212.93 


2213.58 


81 


25 


< 13.20 




2 366 




9.0 


1550.77 


2216.61 


2217.26 


< 50 




< 13.44 








-1.6 


1215.67 


1737.63 


1738.14 


< 269 




< 13.97 








2.0 


977.02 


1396.51 


1396.92 


< 259 




< 13.90" 






0.47436 




1548.19 


2281.82 


2283.46 


645 


37 


> 14.55 




1 511 




-6.2 


1550.77 


2285.62 


2287.25 


524 


33 


> 14.75 








105.2 


1215.67 


1790.54 


1795.92 


2768 


327 


> 14.71* 








7.3 


1206.50 


1778.44 


1779.41 


581 


126 


> 13.43* 








8.1 


1260.42 


1857.67 


1859.04 


617 


31 


> 14.03 








-0.0 


1393.76 


2054.25 


2055.56 


635 


52 


> 13.84* 








-1.7 


1402.77 


2067.54 


2068.85 


475 


56 


> 14.02* 






0.48328 




1548.19 


2295.80 


2296.92 


246 


27 


14.04 


0.08 


1 511 




-3.4 


1550.77 


2299.61 


2300.74 


183 


29 


14.12 


0.09 






28.5 


1215.67 


1802.80 


1803.92 


599 


105 


> 14.04" 








-0.2 


1206.50 


1789.01 


1790.12 


631 


124 


> 13.47* 








0.4 


1260.42 


1869.08 


1869.97 


233 


21 


> 13.52 







HE0515-4414 (zem = 1.710) 

0.50601 ~ 1548.19 2330.91 2332.37 448 28 > 14.36 ~ 2 482 

-0.8 1550.77 2334.79 2335.94 333 21 14.39 0.04 
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Table 4 — Continued 



(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) 

21548 "J'fots <^Wr log N (Tlog N G Flag 

(kms-1) (A) (A) (A) (mA) (mA) 



0.73082 • • • 1548.19 2679.29 2679.99 27 9 12.86 0.14 2 354 

-0.1 1550.77 2683.74 2684.44 < 17 ■•■ < 12.96 

0.94042 ••• 1548.19 3003.24 3005.31 365 20 > 14.27 ■•• 1 499 

2.1 1550.77 3008.24 3010.31 242 23 14.24 0.05 

32.9 1215.67 2357.83 2360.83 1052 29 > 14.68 

20.0 1206.50 2340.48 2342.09 190 18 13.08 0.05 

10.9 1260.42 2445.09 2446.77 122 19 13.00 0.07 

6.9 1393.76 2703.77 2705.23 233 13 13.58 0.03 

7.3 1402.77 2721.26 2722.73 56 11 13.15 0.08 



HS0624+6907 (zem = 0.370) 

0.06351 ~ 1548.19 1646.18 1646.78 106 10 Wji 004 1 503 

-7.7 1550.77 1648.92 1649.52 61 10 13.55 0.07 

-16.2 1215.67 1292.28 1293.28 594 11 > 14.43 

-20.7 1025.72 1090.36 1091.20 357 23 > 15.09 

-1.6 1206.50 1282.91 1283.30 150 9 13.05 0.04 

0.07574 ■•• 1548.19 1665.24 1665.72 53 6 13.19 0.05 1 503 

5.4 1550.77 1668.01 1668.49 32 8 13.23 0.11 

0.9 1215.67 1307.45 1308.15 300 10 > 14.09 

1.3 1025.72 1103.16 1103.68 127 28 14.49 0.12 



HS0747+4259 (ze„ = 1.897) 

0.83662 ~ 1548.19 2842.55 2844.22 303 25 > 14.11 ~ 2 481 

-0.7 1550.77 2847.27 2848.95 232 30 14.20 0.06 

-19.8 1215.67 2232.02 2233.33 < 86 • • ■ < 13.39 



HS0810+2554 (zem = 1.510) 

0.83135 ~ 1548.19 2834.27 2836.40 774 49 > 14.63 ~ 1 499 

2.8 1550.77 2838.99 2841.12 724 53 > 14.89 

62.4 1215.67 2225.35 2228.26 1242 106 > 14.61 

36.6 1206.50 2209.27 2210.36 380 62 > 13.25' 

13.1 1393.76 2551.82 2553.29 262 30 13.67 0.08 

13.0 1402.77 2568.32 2569.80 271 31 13.91 0.06 

0.87687 ••• 1548.19 2905.00 2906.32 219 55 > 13.73' ■■• 2 486 

-7.8 1550.77 2909.83 2911.15 162 49 < 13.94 

-24.5 1215.67 2281.20 2282.09 < 98 • • • < 0.30 



PG0953-F415 (zem = 0.234) 

0.06807 ~ 1548.19 1653.40 1653.77 136 26 > 13.53» ~ 1 383 

-2.7 1550.77 1656.15 1656.52 < 50 •■• < 13.63 

3.7 1215.67 1298.15 1298.74 281 8 > 14.08 

-12.3 1025.72 1095.32 1095.81 129 8 14.44 0.03 

-0.2 977.02 1043.41 1043.65 74 7 13.21 0.05 

-3.1 1031.93 1102.00 1102.31 108 7 14.13 0.03 

-5.1 1037.62 1108.08 1108.39 74 7 14.27 0.04 

-3.4 1238.82 1322.98 1323.29 48 7 13.49 0.06 

3.0 1242.80 1327.24 1327.5.-; 29 7 13.49 0.10 



MARK132 (zem = 1.757) 

0.70776 ~ 1548.19 2643.55 2644.28 29 9 TTSS 013 2 354 

0.4 1550.77 2647.95 2648.68 < 17 ■ • ■ < 12.93 

0.74843 ••• 1548.19 2706.10 2707.69 318 13 14.13 0.03 1 482 

-0.4 1550.77 2710.60 2712.20 165 12 14.08 0.03 

0.76352 ••• 1548.19 2729.63 2730.92 110 11 13.53 0.05 1 482 

-1.4 1550.77 2734.17 2735.46 75 16 13.64 0.09 
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(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


^1548 






A, 


Ah 






log AT 


flog N 


G 


Flag 




(kms-1) 


(A) 


(A) 


(A) 


(mA) 


(mA) 




















3C249.1 (zo 


„ = 0.312) 










0.02616 




1548.19 


1588.55 


1588.81 


24 


7 


12.81 


0.14 


2 


375 




-7.6 


1550.77 


1591.19 


1591.45 


< 13 




< 12.81 










-33.4 


1215.67 


1246.76 


1247.75 


317 


15 


13.95 


0.03 








-31.2 


977.02 


1002.15 


1002.74 


524 


78 


> 13.91 














PG1206-I-459 ( 


em 


1.155) 










0.60072 




1548.19 


2477.60 


2478.89 


185 


25 


13.74 


0.07 


2 


358 




-5.4 


1550.77 


2481.72 


2483.01 


< 45 




< 13.38 








0.73377 




1548.19 


2683.80 


2684.56 


176 


7 


13.80 


0.02 


2 


494 




-2.1 


1550.77 


2688.26 


2689.02 


83 


9 


13.68 


0.05 






0.92677 




1548.19 


2979.54 


2985.83 


2363 


59 


> 15.15 




1 


477 




-75.8 


1550.77 


2984.50 


2990.79 


2156 


54 


> 15.39 










-15.0 


1215.67 


2339.49 


2345.04 


2384 


30 


> 15.07 










66.8 


1206.50 


2321.60 


2326.84 


983 


40 


> 13.97 










-94.8 


1238.82 


2384.15 


2388.91 


876 


39 


> 14.85 










-98.9 


1242.80 


2391.82 


2396.59 


496 


34 


14.82 


0.03 








100.2 


1260.42 


2427.17 


2430.55 


431 


25 


> 13.73 










21.7 


1393.76 


2682.52 


2687.93 


778 


18 


14.08 


0.01 








58.3 


1402.77 


2699.87 


2705.31 


402 


28 


14.06 


0.03 






0.93425 




1548.19 


2993.93 


2995.25 


259 


22 


> 14.18 




1 


511 




0.5 


1550.77 


2998.91 


3000.23 


200 


20 


> 14.34 










-14.4 


1215.67 


2350.48 


2352.09 


552 


21 


> 14.42 










2.6 


1206.50 


2333.27 


2334.11 


133 


15 


> 13.13 










5.9 


1393.76 


2695.43 


2696.40 


143 


11 


> 13.51 










3.6 


1402.77 


2712.87 


2713.84 


98 


11 


13.60 


0.06 












PG1211 + 143 ( 


Zem — 


0.081) 










0.05114 




1548.19 


1626.83 


1627.67 


264 


10 


> 14.07 




1 


509 




3.8 


1550.77 


1629.54 


1630.38 


147 


10 


14.01 


0.03 








36.0 


1215.67 


1276.98 


1279.19 


1214 


7 


> 14.74 










20.8 


1025.72 


1077.59 


1078.89 


823 


103 


> 15.33 










-50.1 


977.02 


1026.58 


1027.04 


304 


85 


> 13.67* 










13.5 


1206.50 


1268.07 


1268.37 


123 


3 


12.94 


0.01 








21.6 


1260.42 


1324.88 


1325.05 


24 


2 


12.29 


0.04 








12.0 


1393.76 


1464.88 


1465.29 


66 


4 


12.96 


0.03 








12.3 


1402.77 


1474.35 


1474.76 


44 


5 


13.04 


0.05 






0.06439 




1548.19 


1647.56 


1648.21 


71 


9 


13.29 


0.05 


1 


511 




-6.9 


1550.77 


1650.30 


1650.95 


31 


9 


13.24 


0.12 








36.4 


1215.67 


1293.48 


1295.26 


861 


5 


> 14.55 










30.6 


1025.72 


1091.36 


1092.50 


448 


12 


> 15.14 










-6.4 


977.02 


1039.69 


1040.14 


144 


9 


13.52 


0.03 








58.3 


1031.93 


1098.13 


1099.07 


183 


11 


14.26 


0.03 








60.6 


1037.62 


1104.19 


1105.13 


67 


11 


14.09 


0.07 








-4.4 


1206.50 


1283.96 


1284.43 


29 


4 


12.17 


0.06 








38.5 


1260.42 


1341.58 


1341.96 


10 


3 


11.85 


0.12 














MRK205 (ze 


™ = 0.071) 










0.00427 




1548.19 


1554.46 


1555.20 


292 


23 


> 14.14 




1 


511 




8.6 


1550.77 


1557.05 


1557.78 


196 


23 


> 14.17 










11.3 


1215.67 


1220.38 


1221.47 


807 


32 


> 14.53 










13.3 


1025.72 


1029.68 


1030.64 


486 


29 


> 15.23 










-6.7 


977.02 


980.81 


981.56 


382 


46 


> 14.13 










4.6 


989.80 


993.68 


994.38 


301 


41 


> 14.51" 










8.4 


1206.50 


1211.42 


1211.97 


263 


32 


> 13.09" 
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(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


21548 


iVabs 




A; 


Ah 






log AT 


flog N 


G 


Flag 




(kms-^) 


(A) 


(A) 


(A) 


(mA) 


(mA) 










ft Q 




1265.57 


1266.25 


240 


1 1 
1± 


\ 1 Q KA 


' ' ' 








3.7 




1399.43 


1400.00 


122 


11 


13.25 


0.04 








1.3 


1 A no TT 


1408.48 


1409.05 


75 


13 


13.33 


0.07 












QSO-123050-(-011522 {ze: 


m = 0.117) 










0.00574 




1548.19 


1556.92 


1557.24 


65 


11 


13.34 


0.07 


1 


511 






1550.77 


1559.51 


1559.75 


36 


Q 

y 


13.33 


nil 








-9.0 


1215.67 


1222.14 


1223.32 


684 


21 


> 14.48 










Q O 

-y.y 


1 onR KC\ 


1213.19 


1213.59 


175 


o4 


^ 1 Q1 ^ 










-12.6 


1 ot?n AO 
iZfDU.4z 


1267.41 


1267.77 


74 


1 1 


12.83 


0.07 








- 1.0 


I'^O.'^ . 70 


1101.17 


1 ini.')2 


Gfi 




12 , 9{) 


0.0) 








-7.9 


1 i(J2. 77 


1111).-")! 


1 IKl.il!) 




i) 


l.'i,21 


().()(> 












PG1241 + 176 ( 


^cm — 


1.283) 










0.48472 




1548.19 


2298.17 


2299.23 


305 


46 


> 13.88" 




2 


352 




-8.9 


1550.77 


2302.00 


2303.05 


< 93 




< 14.15 








0.55070 




1548.19 


2399.44 


2402.21 


852 


44 


> 14.60 




1 


450 




12.8 


1550.77 


2403.43 


2406.21 


580 


50 


14.60 


0.05 






0.55842 




1548.19 


2412.41 


2413.21 


230 


33 


> 13.75'* 




1 


482 




-4.9 


1550.77 


2416.42 


2417.23 


204 


29 


14.26 


0.10 






0.75776 




1548.19 


2721.08 


2721.63 


100 


13 


13.58 


0.07 


1 


486 




-0.1 


1550.77 


2725.61 


2726.06 


50 


12 


13.47 


0.11 






0.78567 




1548.19 


2763.64 


2765.44 


134 


14 


13.58 


0.05 


1 


486 




1 A 


1550.77 


2768.23 


2770.04 


71 


13 


13.59 


0.08 






0.89546 




1548.19 


2934.10 


2934.86 


106 


25 


> 13.42" 




1 


509 




-1 5 


1550 77 


2938.98 


2939.74 


126 


27 


> 13 79" 










-5.6 


1215.67 


2303.47 


2304.94 


510 


42 


> 13.9/ 










-3.4 




2286.53 


2287.17 


79 


23 


< 12.76 














PG1248-I-401 ( 


Zem — 


1.032) 










0.55277 




1548.19 


2403.65 


2404.35 


52 


15 


13.16 


0.13 


2 


486 




8 5 


1550 77 


2407.64 


2408.35 


64 


17 


13 56 


12 






0.56484 




1548.19 


2422.25 


2423.16 


97 


14 


13.51 


0.06 


1 


484 




4.9 


1550.77 


2426.27 


2427.19 


80 


13 


13.71 


0.07 






0.70104 




1548.19 


2633.11 


2633.97 


108 


24 


13.64 


0.13 


1 


492 




-0.3 


1550.77 


2637.49 


2638. S.I 


76 


22 


13.72 


0.13 






0.77291 




1548.19 


2743.91 


2745.76 


619 


33 


> 14.56 




1 


495 




3.7 


1550.77 


2748.48 


2750.33 


564 


29 


> 14.77 










6.2 


1393.76 


2470.19 


2471.84 


426 


14 


> 13.98 










4.7 


1402.77 


2486.17 


2487.83 


315 


12 


14.02 


0.02 






0.85508 




1548.19 


2870.57 


2873.99 


826 


39 


> 14.59 




1 


463 




-11.5 


1550.77 


2875.34 


2878.77 


567 


41 


> 14.65 










-79.1 


1260.42 


2336.91 


2338.58 


196 


29 


13.34 


0.11 








-33.2 


1393.76 


2584.38 


2585.96 


226 


22 


13.52 


0.05 








-45.9 


1402.77 


2601.10 


2602.69 


137 


32 


13.59 


0.10 












PG1259-F593 ( 


Zem — 


0.478) 










0.04615 




1548.19 


1619.46 


1619.85 


103 


21 


> 13.41" 




1 


499 




0.9 


1550.77 


1622.16 


1622.54 


74 


19 


13.71 


0.11 








34.8 


1215.67 


1271.17 


1272.78 


1043 


18 


> 14.65 










-11.3 


1025.72 


1072.39 


1073.64 


748 


35 


> 15.40 










-7.2 


977.02 


1021.61 


1022.67 


411 


22 


> 14.14 










35.0 


1206.50 


1261.91 


1262.83 


220 


13 


13.14 


0.03 







PKS1302-102 (zem = 0.278) 
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Table 4 Continued 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) (11) 


21548 




A,, 


A, 


A„ 


W,. 




log AT 


flog N 


G Flag 




(kms^^) 


(A) 


(A) 


(A) 


(niA) 


(mA) 










0.00438 




1548.19 


1554.81 


1555.14 


32 


8 


12.94 


0.11 


1 383 




-6.5 


1550.77 


1557.40 


1557.73 


< 18 




< 12.98 








-1.3 


1215.67 


1220.75 


1221.23 


300 


18 


> 14.15 








-0.1 


1025.72 


1029.83 


1030.58 


240 


14 


> 14.88 








-17.3 


977.02 


980.89 


981.60 


234 


34 


> 13.91 






CS0873 (Zem = 1.014) 


0.66089 




1548.19 


2569.87 


2572.23 


424 


37 


> 14.24 




1 494 




3.7 


1550.77 


2574.14 


2576.51 


272 


34 


> 14.31 






0.73385 




1548.19 


2683.86 


2684.74 


77 


19 


13.38 


0.10 


1 366 




-2.2 


1550.77 


2688.33 


2689.21 


< 39 




< 13.33 












PG1630+377 ( 




1 476) 








0.75420 




1548.19 


2715.42 


2716.22 


49 


14 


13.17 


0.11 


2 354 




-1.2 


1550.77 


2719.94 


2720.74 


< 27 




< 13.15 






0.91440 




1548.19 


2963.34 


2965.06 


384 


16 


> 14.28 




1 469 




10.3 


1550.77 


2968.17 


2969.99 


317 


18 


> 14.43 








-2.0 


1215.67 


2325.82 


2328.60 


1076 


33 


> 14.69 








8.9 


1206.50 


2309.25 


2310.59 


142 


19 


12.99 


0.07 






-12.6 


1393.76 


2667.64 


2668.60 


123 


15 


13.26 


0.06 






-12.5 


1402.77 


2684.89 


2685.86 


98 


12 


13.47 


0.06 




0.95269 




1548.19 


3022.00 


3023.84 


560 


21 


> 14.43 




1 477 




-24.3 


1550.77 


3027.02 


3028.87 


276 


22 


14.27 


0.04 






70.6 


1215.67 


2372.75 


2376.16 


1434 


25 


> 14.83 








8.9 


1206.50 


2355.41 


2356.56 


157 


19 


13.00 


0.06 






3.9 


1260.42 


2460.75 


2461.78 


43 


13 


12.54 


0.13 






14.7 


1393.76 


2720.92 


2722.40 


202 


15 


13.49 


0.04 






23.1 


1402.77 


2738.52 


2740.01 


136 


21 


13.59 


0.07 










PG1634+706 ( 




1.334) 








0.41935 




1548.19 


2197.10 


2197.71 


- 


11 


13.07 


0.11 


2 354 




0.4 


1550.77 


2200.76 


2201.37 


< 22 




< 13.06 






0.65126 




1548.19 


2556.26 


2556.77 


67 


6 


13.33 


0.04 


1 496 




-3.0 


1550.77 


2560.51 


2560.94 


67 


6 


13.61 


0.04 






0.3 


1215.67 


2007.05 


2007.67 


78 


21 


13.25 


0.12 




0.65355 




1548.19 


2559.31 


2560.62 


393 


9 


> 14.29 




1 499 




2.0 


1550.77 


2563.57 


2564.88 


284 


10 


14.34 


0.02 






0.9 


1215.67 


2009.41 


2010.92 


630 


30 


> 14.45 








-63.3 


1206.50 


1993.91 


1995.48 


502 


28 


> 13.71 








-7.2 


1393.76 


2303.95 


2305.30 


145 


11 


13.30 


0.04 






0.1 


1402.77 


2318.85 


2320.21 


73 


10 


13.26 


0.06 




0.81814 




1548.19 


2814.55 


2815.30 


90 


4 


13.51 


0.02 


1 501 




-0.1 


1550.77 


2819.23 


2819.98 


67 


4 


13.61 


0.03 






-0.2 


1215.67 


2209.74 


2210.75 


201 


10 


> 13.92 








-23.7 


1206.50 


2193.22 


2193.83 


115 


7 


12.88 


0.03 






-2.4 


1260.42 


2291.25 


2291.88 


21 


6 


12.21 


0.13 






-1.0 


1393.76 


2533.54 


2534.40 


< 15 




< 12.25 








-3.4 


1402.77 


2549.93 


2550.79 


21 


7 


12.71 


0.13 




0.90560 




1548.19 


2949.76 


2950.74 


198 


7 


14.04 


0.04 


1 509 




-0.5 


1550.77 


2954.66 


2955.65 


165 


7 


14.20 


0.03 






-1.6 


1215.67 


2316.12 


2317.05 


315 


6 


> 14.16 








-5.3 


1025.72 


1954.23 


1955.01 


180 


20 


> 14.39* 








-4.3 


1031.93 


1966.01 


1966.81 


104 


22 


14.08 


0.10 






-7.1 


1037.62 


1976.85 


1977.66 


67 


18 


14.11 


0.12 






-2.7 


1206.50 


2298.85 


2299.38 


61 


6 


12.58 


0.05 






1.3 


1238.82 


2360.43 


2360.98 


67 


6 


13.57 


0.04 
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Table 4 — Continued 





^9^ 
' 




(4) 


(5) 




('7\ 


\°) 




(iOJ 




^1548 


^Vabs 




Ai 


Ah 


T/T/ 
yv r 




log iV 




c 


Flag 




(km s ) 


(A) 


(A) 


(A) 


(inA) 


(inA) 












-8.8 


1242.80 


2368.02 


2368.57 


52 


5 


13.74 


0.05 








-21.3 


1260.42 


2401.59 


2402.15 


107 


4 


13.09 


0.03 








7.0 


1393.76 


2655.64 


2656.26 


259 


5 


13.82 


0.03 








-3.1 


1402.77 


2672.82 


2673.44 


40 


6 


13.01 


0.06 






0.91144 




1548.19 


2958.77 


2959.74 


34 


7 


12.96 


0.09 


2 


383 




-2.8 


1550.77 


2963.69 


2964.66 


< 17 




< 12.92 










18.3 


1215.67 


2322.48 


2324.84 


660 


11 


> 14.48 










1.6 


1025.72 


1959.76 


1961.33 


311 


25 


> 14.91 














HS1700+6416 


Zem — 


2.736) 










0.08077 




1548.19 


1673.11 


1673.40 


289 


90 


> 13.80 




2 


371 




-1.5 


1550.77 


1675.90 


1676.18 


< 154 




< 14.40 










4.7 


1215.67 


1313.28 


1314.36 


661 


28 


> 14.43 














PG1718+481 ( 


^em — 


1.084) 














1548 19 


2259.21 


2259.97 


147 


13 


13 68 


04 


1 


480 




-8 6 


1550 77 


2262.97 


2263.72 


120 


13 


13 83 


05 












H1821+643 


ETTl — 


0.297) 










0.22503 




1548.19 


1895.74 


1897.42 


137 


23 


13.61 


0.07 


1 


381 




3.0 


1550.77 


1898.90 


1900.57 


< 47 




< 13.38 










-51.7 


1215.67 


1488.23 


1489.74 


834 


13 


> 14.60 










-33.7 


1025.72 


1255.90 


1256.80 


496 


9 


> 15.24 










-17.9 


977.02 


1196.50 


1197.22 


306 


13 


> 13.99 










-11.5 


1031.93 


1263.72 


1264.45 


172 


9 


14.26 


0.02 








-12.0 


1037.62 


1270.68 


1271.42 


106 


9 


14.30 


0.04 








-4.0 


1206.50 


1477.75 


1478.19 


88 


7 


12.79 


0.04 






0.24531 




1548.19 


1927.27 


1928.76 


75 


22 


13.32 


0.12 


2 


383 




-1.2 


1550.77 


1930.47 


1931.97 


< 44 




< 13.35 










3.2 


1215.67 


1513.68 


1514.08 


47 


8 


13.01 


0.07 








-25.3 


1025.72 


1277.01 


1277.51 


55 


8 


13.95 


0.06 








1.1 


1031.93 


1284.86 


1285.28 


50 


7 


13.67 


0.06 








3.8 


1037.62 


1291.94 


1292.37 


29 


5 


13.70 


0.07 














PHL1811 (zc 


,„ = 0.190) 










0.08091 




1548.19 


1673.01 


1673.74 


189 


25 


> 13.93 




1 


511 




-0.6 


1550.77 


1675.79 


1676.53 


129 


22 


14.03 


0.10 








-27.8 


1215.67 


1313.15 


1314.77 


904 


9 


> 14.64 










-36.0 


1025.72 


1108.19 


1109.04 


622 


8 


> 15.36 










-19.7 


977.02 


1055.61 


1056.31 


338 


9 


> 14.07 










11.7 


989.80 


1069.72 


1070.02 


137 


7 


> 14.47 










10.9 


1206.50 


1303.98 


1304.27 


179 


4 


> 13.31 










2.2 


1260.42 


1362.31 


1362.60 


172 


4 


> 13.44 









^COG 

Note. — C IV systems by target and redshift of C IV 1548. Upper limits are 2cr limits for both and log N . The 
column densities were measured by the AODM, unless "COG" indicated, in which case, the limit is from assuming 
Wj. results from the linear portion of the COG. The definite CIV doublets are labeled group G — 1, while the likely 
doublets are G — 2. The binary flag is described in Table[3]and ^ 13.21 
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3.3. Final C IV Sample 

The 319 candidate C IV systems with WrS548 > 
3aWr, 1548, Rw within the expected range, and 
l^'i'Civl ^ lOkms"^ were visuahy evaluated by 
multiple authors. The remaining systems were 
reviewed by at least one author. We assessed the 
likehhood of each candidate, utihzing the diag- 
nostics described in the previous section. Then 
we agreed upon two groups of intervening C IV 
systems that had 21548 more than 1000 km s^^ 
redward of the Milky Way Galaxy (z — 0) and 
more than 3000 km s~^ blueward of the back- 
ground source {z = Zem)- 44 definite systems 
(G = 1) and 19 likely systems (G = 2). The de- 
tailed information about all G = 1-1-2 systems is 
listed in Table |4] and the summary of the doublets 
is listed in Table [3 An example velocity plot of 
a C IV doublet is given in Figure [H and the full 
complement of velocity plots are in Appendix [A] 
The cumulative distribution of column densities 
and equivalent widths are given in Figure U 

The final, visual evaluation of the candidates 
was subjective. We required consensus among the 
authors in defining our final sample. Typically, the 
definite (G = 1) systems show other transitions, 
usually Lya and/or Si IV. When there was no cov- 
erage of the other transitions, which had lower rest 
wavelengths, the G = 1 C IV doublets were multi- 
component with matching AOD profiles. 

Typically, the G = 2 systems include the C IV 
doublets detected in regions where the S/N was 
low, so the comparison of the profiles was less con- 
clusive, and there were not enough other positive 
diagnostics {e.g., detection of Lya) to boost the 
confidence of the identification. 

For all C IV systems, we confirmed that the 
C IV doublet was not a higher-redshift O I 1302, 
Si II 1304 pair, which has rest wavelength sepa- 
ration similar to that of the doublet, by checking 
whether Si II 1260 existed and was much stronger 
than the 1304 line. We also confirmed that the 
doublet and the other transitions were not other, 
common transitions at different redshifts. 

We adjusted the wavelength Hmits (A/,Aft) as 
necessary to correct for blending. For example, in 
Figure [2I the wavelength bounds for the C IV 1548 
line was reduced. 

For the remaining analyses, we restrict our in- 
tergalactic C IV sample to the systems where both 



lines of the doublet were detected with Wr > 
ScTw^H This left 38 G = 1 systems and five G 
= 2 systems. 




13.0 13.5 14.0 14.5 15.0 15.5 



3 
o 
Q 



t 
3 

'Z, 



o 
H 




'r,1548 



(mA) 



Fig. 4. — Cumulative column density and equiva- 
lent width distributions for the Wr > Saw^ G = 
1 and G = 2 samples. The black and gray curves 
are the G = 1 and G = 1-1-2 groups, respectively. 
For each absorber, log V(C+'^) and Wr, 154s are in- 
dicated with the hashes across the top, with the 
color indicating the sample. As in Figure [TJ the 
saturated doublets are shown in the middle row 
of hashes in the top panel. The 21548 = 0.92677 
C IV doublet associated with the DLA in the 
PG1206-h459 sightline, with H/^4548 = 2364mA, 
is not shown in the W^r,i548 plot (bottom panel). 



''Wc list all G = 1+2 C IV systems in Table|4]for the com- 
munity and any future analyses with different criteria. 
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Table 5 
C IV DOUBLET SUMMARY 



Target 


G 


(■i) 

21548 


i A\ 
VKr-,1548 

(mA) 


(mA) 


(6) 

log A''l548 




(T) 
log A^'isso 


logAf^C ) 


pnm 1 7-1-91 


2 


u . o -L 


94 ± 24 


< 45 


1 ^ 

lo. 


/1 7 -L n 11 

. 4 / i . 1 1 




1 

lo. 


,oy 


1 Q 47 _|_ n 1 1 




]^ 


n 57R'^9 




442 ± 25 


14. 


56 




14 


52 ± 03 


1 4 CO _|_ n no 


PKS0232— 04 


1 


0.73910 


372 =b 34 


359 ±41 


14 


. ±D ^ U . Uo 




14 


43 


14. 1 c + f) no 




2 




69 di 22 


< 37 


<, lo. 


13 




1 

lo. 


■^9 
,OZ 


< 13 13 






'^fi071 


86 it 17 


< 35 


1 

lo. 


■^7 -L n nQ 




1 

lo. 


97 
Z ( 


1 Q '^7 -1- n HQ 


PKS0454-22 


2 




142 ± 24 


81 ±26 


1 
lo. 


ft/1 -L n na 




1 

lO, 


7n -L n 1 /I 

. f U ^ U. 14 


1 ftc _|_ n n7 






0.24010 


644 ± 57 


495 ± 54 


> 14, 


ona 
,zu 




14 


.oy 


> 14.39 






0.27797 


274 ± 52 


< 104 


^ lO. 


83^ 




14 


20 


n 3 83 1 4 1 41 






U.OCSiOZ 


Oft -1- QO 

yo It oz 


1 1 n J- 07 

iiU It Z / 


< 13. 


.55 


> 


13, 


,73^ 


n Q y1 Q 10 771 1> 
[10.40, lo. 1 1 \ 




1 


Xj.^KjZZ ( 


1 / f ± zy 


140 ± 31 


< 13. 


,72 


> 


13. 


,84^ 


Fl /I 10 QOl 

[13. o4, I0.S8I 




2 


c\ A on c n 


81 i 25 


< 50 


< 13. 


,20 


< 


13. 


,44 


< 13.44 




1 


n A '74 K 


645 i 37 


524 ± 33 


> 14. 


,55 


> 


14. 


,75 


> 14.75 






n /I S'?9S 

U .4:OOZO 


Z4D i z ; 


1 S"^ -1- 9Q 

loo i zy 


14. 


n A _L n no 

.04 ± U.U8 




14. 


1 o _i_ n nn 

,lz ± 0.U9 


1 A n7 -L n nfi 

14. U 1 ZIZ U.UD 


rlEjUO i O— ^■y: 141: 


2 


U.oUDUi 


A AH -\- 9S 
440 ± ZO 


Q -1- 91 
OoO It Zi 


> 14, 


,36 




14, 


o n _L n n A 

,39 ± 0.04 


1 A QQ -L n nA 

14. oy It U.U4 




2 


U. * oUoz 


27 i 9 


< 17 


12, 


.86 ± 0.14 


< 


12, 


.96 


1 Qft _L n 1 /I 
Iz.oD ih U.14 




1 


U.y4U4Z 


Qf3 K -L on 
oDO ih zU 


o/i o _i_ oo 
Z4Z it Zo 


> 14, 


.27 




14, 


.24 ± 0.05 


1 A OA -L n nK 
14. Z4 it U.Ut) 


rlk3UD^4-|-0yU 1 


± 


U.UDooi 


1 nft -Lin 
iUD It lU 


fti -Lin 
Oi It iU 


13, 


e A 1 A 

.54 ± 0.04 




13, 


.55 ± 0.07 


1 KA -L n no 
1O.04 It U.Uo 




1 


U.U / D / 4 


K'i -4- R 
Do It D 


QO -L fi 
OZ It O 


13, 


.19 ± 0.05 




13, 


.23 ± 0.11 


1 Q on -L n CiA 

lO.ZU It U.U4 


u Q n 7 /1 7 _i_ /I o Q 
rioU / 4 / -|-4zoy 


2 


U.OODDZ 


T m -1- 9 r 
oUo It 


9^9 -U Qn 
ZoZ It ou 


> 14. 


,11 




14. 


on 1 n nc 

,zO ± 0.U6 


1 A on -L n nft 
14. ZU It U.UD 


rloUoiU+^oo4 


1 


u.ooloo 


( /4 ± 4y 


724 ± 53 


> 14. 


,63 


> 


14. 


,89 


> 14.89 




2 


U.o f DO / 


219 i 55 


ICO _L A n 
162 ih 4y 


> 13. 


,73^ 


< 


13. 


,94 


r 1 Q "7Q 1 A 111 

[13. 10. 14.11J 


TD ("I no _i_ /1 1 
U y O o -|- 4 i 


1 


U.UDoU / 


1 Qfi -1- 9fi 
ioD It ZD 


< 50 


> 13. 


,53'^ 


< 


13. 


,63 


r 1 Q 1 771 
[lo.oo, lo. ( 1 j 


iVi AlxiV ± O Z 


2 


U. / U ( ID 


29 i 9 


< i ( 


12. 




< 


12, 


.93 


1 9 _L n 1 
IZ.OO It U. lo 




1 


U. 1 4o40 


Q1 ft -1- 1 Q 
Oio It lo 


1 Rf; -L 1 O 
IDO It iZ 


14, 


1 o _l_ n no 




14, 


no _L n no 
.08 ih O.Oo 


1/1 1 n -L n no 

14. lU It U.UZ 




1 


U. * DOOZ 


110 i 11 


75 ± 16 


13, 


.53 ± 0.05 




13, 


A In nn 

.64 ± 0.09 


1 Q -L n n/1 
lo. 00 ih U.U4 


Q(^0/1Q 1 

oo^^y. i 




U.UZDlO 


OA -4- T 
Z4 It ( 




12, 


.81 it 0.14 


< 


12, 


.81 


10 (21 -L n 1/1 

IZ.Ol ± U. 14 


Jr lor 1 ZUD y 


n 
Z 


u.ouu / z 


1 ftp; -1- op; 

loo It ZD 


^ AK 
<. 40 


13, 


.74 ± U.U7 


< 


13, 


.38 


1 q 7/1 _L n n7 
lo. ( 4 It U.U 1 




o 


U. ( oo 1 / 


i / D It / 


83 ± 9 


13. 


.oU ± U.Uz 




13. 


.d8 ± 0.05 


1 77 _L n no 
lo. / 1 It U.UZ 




1 


U.yzD f / 


zoDo ± oy 


zlOD ih o4 


> 15. 


.15 


> 


15. 


.39 


> 15.39 




1 


n HQ /I n c; 
U.yo4zo 


O _L OO 

zoy ± zz 


onn _L on 
ZUU it ZU 


> 14. 


.18 


> 


14. 


.34 


> 14.34 


■pr^ 1 01 1 1 1 /I Q 


± 


n n ^ 1 1 /I 


ofi/i -Lin 
zd4 It iU 


147 ± 10 


> 14. 


07 




14. 


.01 ±0.03 


1/1 m -Ln m 

14. Ul ± U.Uo 




1 
1 


U.UD4oy 


/ i It y 


31 ± 9 


13. 


.29 ±0.05 




13. 


.24 ±0.12 


1 Q 9s -L n nr^ 
lo.Zo It U.Uo 




1 


U.UU4Z ( 


zyz ± Zo 


1 Qd -L OQ 

lyo It Zo 


> 14, 


.14 


> 


14, 


.17 


\ 1 /I 17 
14, 1 ( 


r^Q/^ 1 OOHKH 1 m 1 KOO 

Ci^oU— izoUoU+Ui iozz 


1 


U.UUo /4 


65 i 11 


36 ± 9 


13, 


.34 ±0.07 




13, 


.33±0.11 


1 Q Oyi -L n nft 
lo.o4 lb U.Uo 


PG1241+176 


2 


0.48472 


305 ± 46 


< 93 


> 13, 


.88'' 


< 


14, 


.15 


[13.88,14.01] 




1 


0.55070 


852 ± 44 


580 ± 50 


> 14, 


.60 




14, 


.60 ±0.05 


14.60 ±0.05 




1 


0.55842 


230 ± 33 


204 ± 29 


> 13. 


.75^ 




14. 


.26 ±0.10 


14.26 ±0.10 




1 


0.75776 


100 ± 13 


50± 12 


13. 


.58 ±0.07 




13. 


.47 ±0.11 


13.54 ±0.06 




1 


0.78567 


134 ± 14 


71 ± 13 


13. 


.58 ±0.05 




13. 


.59 ±0.08 


13.58 ±0.04 




1 


0.89546 


106 ±25 


126 ± 27 


> 13. 


.42° 


> 


13. 


.79° 


> 13.79 


PG1248+401 


2 


0.55277 


52± 15 


64 ±17 


13. 


.16±0.13 




13. 


.56 ±0.12 


13.25 ±0.09 




1 


0.56484 


97 ±14 


80 ±13 


13, 


.51 ±0.06 




13, 


.71 ±0.07 


13.57 ±0.05 




1 


0.70104 


108 ± 24 


76 ±22 


13, 


.64 ±0.13 




13, 


.72 ±0.13 


13.68 ±0.09 




1 


0.77291 


619 ±33 


564 ± 29 


> 14, 


.56 


> 


14. 


.77 


> 14.77 




1 


0.85508 


826 ± 39 


567 ± 41 


> 14. 


.59 


> 


14. 


.65 


> 14.65 


PG1259+593 


1 


0.04615 


103 ±21 


74 ±19 


> 13. 


.41° 




13. 


.71 ±0.11 


13.71±0.11 


PKS1302-102 


1 


0.00438 


32 ±8 


< 18 


12. 


.94 ±0.11 


< 


12. 


.98 


12.94±0.11 


CS0873 


1 


0.66089 


424 ±37 


272 ± 34 


> 14. 


.24 


> 


14. 


.31 


> 14.31 




1 


0.73385 


77± 19 


< 39 


13. 


.38 ± 0.10 


< 


13. 


.33 


13.38 ±0.10 


PG1630+377 


2 


0.75420 


49± 14 


< 27 


13. 


.17±0.11 


< 


13. 


.15 


13.17±0.11 




1 


0.91440 


384 ± 16 


317 ±18 


> 14, 


.28 


> 


14, 


.43 


> 14.43 




1 


0.95269 


560 ±21 


276 ± 22 


> 14, 


.43 




14, 


.27 ±0.04 


14.27 ±0.04 


PG1634+706 


2 


0.41935 


44± 11 


< 22 


13. 


.07±0.11 


< 


13. 


,06 


13.07±0.11 




1 


0.65126 


67 ±G 


67 ±6 


13. 


.33 ±0.04 




13. 


.61 ±0.04 


13.40 ±0.03 




1 


0.65355 


393 ±9 


284 ± 10 


> 14, 


.29 




14, 


.34 ±0.02 


14.34 ±0.02 




1 


0.81814 


90 ±4 


67 ±4 


13, 


.51 ±0.02 




13, 


.61 ±0.03 


13.54 ±0.02 
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4. Survey Sensitivity 

4.1. Monte-Carlo Completeness Tests 

We used Monte-Carlo tests to determine the 
column density and equivalent width limit, as a 
function of the redshift, where we recovered 95% 
of the simulated C IV doublets for each spec- 
trum. From these sensitivity functions, we esti- 
mated the unblocked co-moving absorption path- 
lengths AX{N{C+^)) and AX(VKr,i548) for our 
survey, as defined below. 

The simulated C IV doublets were single- or 
multi-component Voigt profiles. Each Voigt pro- 
file was defined by a randomly selected redshift, 
column density, and Doppler parameter. For 
the multi-component simulated doublets, each 
component was assigned a decreasing fraction 
of the total column density and a velocity off- 
set from the assigned redshift. For each col- 
umn density bin, AlogA^(C+^) — 0.1, in the 
range 12.7 < logN{C+^) < 14.3, one-hundred 
randomly-generated Voigt profiles were added to 
the spectrum for each redshift bin, Sz — 0.005, 
in the range covered by the spectrum. The one- 
hundred simulated doublets were distributed in 
a manner that avoided blending. The simulated 
doublets had Doppler parameters 5kms~^ < 6 < 
25kms^^ and from one to five components in the 
range — lOOkms""'^ < Svabs < 100 kms^^. We vi- 
sually inspected the synthetic doublets to verify 
our model profiles resembled observations of z > 2 
C IV doublets. 

The profiles were added to the "cleaned" spec- 
trum with the appropriate noise. The cleaned 
spectrum was the original spectrum with all 
automatically-detected absorption features re- 
placed with continuum and noise. The noise was 
randomly sampled from the neighboring pixels, 
not in the absorption features. We elected to mod- 
ify the original spectra (as opposed to generating 
completely synthetic spectra) to keep the pixel- 
by-pixel signal-to-noise properties of the spectra 
realistic. The minimum fiux level (i.e., the fiux 
at line black) for the optically-thick Voigt profiles 
was measured empirically from the troughs of the 
strongest absorption lines, typically Galactic, in 
the original spectrum. 

The blind doublet search described in § 13. II was 
conducted on the synthetic spectra. The recovered 



doublets were matched to the input simulated dou- 
blets. The 95% completeness limits log Niijn{C~^^) 
and VFr,iim (for the 1548 line) were measured for 
each redshift bin. For each redshift bin, we fit 
a fourth-order polynomial to the percentage of 
simulated C IV doublets recovered and solved for 
logA^hm(C+^) and W^rjim (see Figures [5] and El 
respectively). 

4.2. Search Pathlength 

The sensitivity functions were constructed for 
every sightline by the previously described Monte- 
Carlo method. From these, we measured the co- 
moving absorption pathlength AX{N{C~^^)) and 
AX(Wr,i548) (see Figure We converted the 
sensitivity functions from redshift space {Sz = 
0.005 from the Monte-Carlo tests) to pathlength 
space as follows: 

x{z) = --^^nM{i + z)3 + n^ (9) 

6X{z) = X{z + 0.5Sz)-X{z-0.5Sz), 

where ilM = 0-3 and Qjy — 0.7 for our adopted 
cosmology. The total pathlength as a function of 
column density AX{N{C^^)) is the sum of the 
bins SX where the 95% completeness limit is < 
iV(C+3): 

AX{N{G+^)) = '5^(N^) (10) 

Nii„<Ni 

The error a ax refiects the uncertainty induced by 
the uncertainty in A^(C+^) and the interpolation of 
due to the column density bins i5N. The afore- 
mentioned also applies to the pathlength as a func- 
tion of equivalent width AX(VFr,i548)- For exam- 
ple, in Figure [H the bins where logiVHm(C+'^) < 
13.5 totaled AX{N{C+^)) = 0.44 for CS0873 
(zem — 1.014). Similarly, we found AX(Wr, 1548) = 
0.44 for Wr.vnu < 83mA. The total absorption 
pathlength, 'AX(iV(C+3)) (or AX(W^r,i548)) is 
the sum (Equation [TU)) over all sightlines. 

We masked out common Galactic lines, which 
we determined from stacked STIS E140M, E230M, 
and G230M spectra. We measured the wavelength 
bounds (XiiXh) of the common Galactic lines in 
the stacked spectra. We buffered the redshift 
range excluded due to Galactic hues by ±20% of 
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Table 5 — Continued 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 




(7) 


(8) 


Target 


G 


^1548 


Wr,1548 




log TV 1548 




log Afi550 


logTV(C+^) 








(mA) 


(mA) 












1 


0.90560 


198 ± 7 


165 ±7 


14.04 ±0.04 




14.20 ±0.03 


14.12 ± 0.03 




2 


0.91144 


34 ± 7 


< 17 


12.96±0.09 


< 


12.92 


12.96±0.09 


HS1700+6416 


2 


0.08077 


289 ± 90 


< 154 


> 13.80 


< 


14.40 


[13.80, 14.06] 


PG1718+481 


1 


0.45953 


147 ± 13 


120 ±13 


13.68 ±0.04 




13.83 ±0.05 


13.73 ±0.03 


H1821+643 


1 


0.22503 


137 ± 23 


< 47 


13.61 ±0.07 


< 


13.38 


13.61 ±0.07 




2 


0.24531 


75 ±22 


< 44 


13.32 ± 0.12 


< 


13.35 


13.32 ± 0.12 


PHL1811 


1 


0.08091 


189 ± 25 


129 ± 22 


> 13.93 




14.03 ±0.10 


14.03 ±0.10 



'^log measured by assuming Wr results from the linear portion of the COG. 
^imits set by doublet lines increased/decreased by la to set log A(C"'"^) range. 



Note. — Summary of CIV doublets by target and redshift of CIV 1548. The definite CIV doublets are labeled 
group G — 1, while the likely doublets are G — 2. Upper limits are 2cr limits for both Wr and log A. The adopted 
column density for the CIV doublets are listed in the last column (see 12. 4t . 



the width of the hnes (see Figures [S] and O ■ Each 
Galactic hne affected the search path twice, as 
a contaminant to both C IV hnes. We excluded 
the regions 1000 km s^^ redward of the Galaxy 
and 3000 km s^^ blueward of the the background 
source. 

We also masked out the saturated regions in 
the Lya forest for each spectrum. We excluded 
pixels with f\i < Uf^i, and ±3 neighboring pixels, 
in all features detected by the automatic feature- 
finding algorithm (see ij 13. ip . The excluded pixels 
included ones from the strongest C IV doublets, 
but this amounted to a small fraction (< 1%) of 
the total AX. 

4.3. Lya Forest Contamination 

The Lya forest was the largest contamination 
in the C IV doublet survey. Approximately 60% of 
the total pathlength resides in the Lya forest. For 
most of these systems, the detection of associated 
transitions such as Lya absorption or the Si IV 
doublet lent credibility to the C IV identification. 
For several of our systems, however, we did not 
have coverage of other transitions. An example 
of this scenario is the G = 1 doublet at Z1548 = 
0.74843 in the MARK132 sightline (which also had 
a nondescript profile; see Appendix [A}. For cases 
such as this, one must consider the possibility that 
a pair of Lya lines mimicked a C IV doublet by 
chance. Instead of excluding systems where we did 
not have the wavelength coverage or data quality 
to detect associated transitions, we accounted for 



the contamination by considering the rate at which 
Lya transitions masquerade as C IV. 

In principle, the contamination rate could be 
addressed in a variety of ways. While we could 
have considered random lines of sight through cos- 
mological simulations, our purpose was first to 
identify if contamination by Lya was significant. 
We thus adopted the simpler approach of gener- 
ating synthetic spectra from the known proper- 
ties of Lya statistics and searched this pure Lya 
spectrum for putative C IV doublets. This ap- 
proach had the advantages that artificial spectra 
were generated very quickly and we could assess 
our contamination in a Monte-Carlo sense. Since 
the higher redshift E230M data were far more sus- 
ceptible to contamination (the Lya forest lies in 
the range 0.8 ;$ ^ ;$ 1.5 in these data), we took 
the E230M wavelength range as our template. 
In order to generate the spec tra, we took the 



statist ics of the Lya forest from iJanknecht et al 
(|200fi[ ). 'n the redshift range of interest, there was 



only very weak evolution of the Lya line density, 
so for simplicity, we adopted a fixed line dcnsitjo 
(Wl^yjdz = 150 for 13 < logA^ni < 16.5. The 
column density and Doppler parameter distribu- 
tions were taken to be a power law and trun- 
cated Ga ussian, respectively, w ith the parame- 



ters from ' Jankn echt et al. ( 2006f l. The placement 
of the Lya lines was somewhat less clear. At 
high redshift, strong Lya absorbers, identified in 



Above z ~ 1.5, the Lya line density evolves strongly (c./. 
iKim et a"Lll200lD 
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Fig. 5. — Example Monte-Carlo completeness test for logNumiC^^). The column density limit where our 
survey was 95% complete is shown for the E230M spectrum of CS0873 (top, left panel). The bottom axis 
is the redshift of the 1548 line, and the top axis is the corresponding co-moving absorption pathlength 
from Equation [9] The general trend follows the S/N profile of the spectrum whose peak sensitivity lies at 
A « 2600 A, corresponding to Z1548 ~ 0.68. The redshift ranges excluded due to Galactic lines are shaded 
(gray); each Galactic line is counted twice, as a contaminant to both C IV lines. The abrupt, narrow spikes 
to infinity are the excluded saturated pixels in the Lya forest. The total pathlength as a function of N{G^^) 
to which the survey is 95% complete, AX{N{C'^^)), is shown for this single sightline (top, right panel). 
For clarity, the gray shaded region lays under the AX {N {C~^^)) curve (also see the total distribution shown 
in Figure [7]). Also shown is the total number of simulated doublets tested (dashed line) compared to the 
number of recovered doublets (solid line) as a function of log 7V(C"'"'^) (bottom panel). 



high-resolution data by the presence of C IV ab- 
sorption, have been shown t o be highly clustered 
( Fernandez-Soto et al ] |l996h . By using heavier 
carbon atoms to trace the underlying gas com- 
ponent structure, those authors showed Lya ab- 
sorption to be highly clustered on small veloc- 
ity scales. In the E230M data, however, we did 



not have the resolution to resolve such small-scale 
structure, and this was reflected ii i the lack of clus- 
tering detected in the UV data of Janknecht et al 



()2006[ ). In order to reproduce the observed trend 



as closely as possible, we thus distributed our Lya 
lines in a uniform fashion. Further, since the ob- 
served statistics did not account for the small-scale 
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Fig. 6. — Example Monte-Carlo completeness test for Wrjim, similar to Figure [5l The rest equivalent width 
limit for the C IV 1548 line where our survey was 95% complete is shown for the E230M spectrum of CS0873 
(top, left panel). The total number of simulated C IV doublets per Wr.isis bin was not constant because 
the simulated doublets were assigned randomly-drawn iV(C^'^) and b (bottom panel). 



blending mentioned above, we imposed a mini- 
mum separation for two lines \Az\ > 7 x 10~^. 
Finally, we added to our spectra the correspond- 
ing higher-order Lyman lines stronger than a nom- 
inal detection level of 15 mA, since they too could 
contribute to the contamination rate. 

The line lists generated in the above fashion, 
which could be generated very quickly, constituted 
the full information available in a "spectrum." We 
thus adopted a Monte-Carlo approach, which was 
appropriate for both low- and high-contamination 
rates. An automated set of criteria was used to 
select pairs of lines that could be misidentified as 
C IV doublets. We selected pairs of lines that lay 



within 30kms~ of the expected position, had 6- 
values within lOkms"^ of each other, and had an 
equivalent- width ratio between 2 : 1 and 1:1. 
These fake doublets were saved for further ex- 
amination, and the simulation was iterated 1000 
times. 

For each misidentified doublet, an actual spec- 
trum was generated, and noise added at the given 
S/N level. All the candidates were visually in- 
spected to exclude clearly mismatched line pairs. 
The final number of accepted doublets was aver- 
aged over the 1000 simulation iterations, yielding 
dA/civ/dX for misidentified doublets. Since lines 
could be quickly generated, and noise easily added. 
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Fig. 7.— Redshift pathlength AX {log N {C+^)) 
and AX(Wr,i548) as a function of C IV 1548 col- 
umn density (top) and equivalent width respec- 
tively (bottom). The black curve is for the full 
redshift {z < 1 sample, the light gray for z < 0.6, 
and the dark gray for z > 0.6. These estimates are 
based on Monte-Carlo analysis and correspond to 
95% completeness limits. 



we assessed the contamination rate at several S /N 
levels that spanned the range of typical values in 
our survey, with the results shown in Figure [51 In 
all cases, the contamination rate was not a signif- 
icant effect, but in § 16. 5[ we discuss several G — 
2 C IV doublets that might be the result of Lya 
forest contamination. 




4 6 8 10 12 14 16 
Signal-to-noise Ratio 



Fig. 8. — Rate of coincident Lya forest lines 
misidentified as C IV doublets as function of S/N. 
The Lya forest at 0.8 < z < 1.5 was simulated in 
a Monte-Carlo fashion, and the incidence of Lya 
forest lines that were "observed" to have similar 
line profiles and wavelength separation as a C IV 
doublet was estimated for a range of S/N. 

5. Analysis 

5.1. Frequency Distributions 

The frequency distribution is the number of ab- 
sorbers Af per column density A^(C+'^) or equiva- 
lent width Wr,i548 bin per absorption pathlength 
sensitive to those absorbers (see § 14. ip : 



AN{C+^)AX{N{C+^)) 



and 



fiWr) 



AA/" 



(12) 



AWrAX{Wr) 

This quantity represents a nearly full description 
of an absorption line survey and is akin to the 
luminosity function as used in galaxy studies. 

We chose to model the frequency distributions 
with power-law functions because it approximated 
the frequency distributions well. The power-law 
functions are of the form: 



/(iV(C+3)) fci4 



/ jV(C+^) 

V ^^0 



+ 3\ \ "N 



and 



f{Wr) = h 



r,0 



(13) 



(14) 
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Table 6 

f{N{C+^)), (Wcw/dX, AND flc^^ SUMMARY 



(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) 

G (z) zi zn M logAr(C + 2) dA/'civ/dz dA/^civ/dX n^+s fci4 an Pks 

(10-**) (cm^) 



Samples by Rcdshift 



1 


0.65355 


0.00574 


0, 


.95269 


36 


(13.15, 16.00) 




4. 


n+1.3 
^-1.0 




0, 


g7+U.l« 

-0.16 


-1 




0.364 




0.65355 


0.00427 


0, 


.95269 


27 


(13.20, 15.39) 




3. 


c + 0.7 

'°-o.e 


> 3.43 












1+2 


0.65355 


0.00574 


0, 


.95269 


41 


(13.15, 16.00) 




4. 


n+1.4 
"-1.1 




0, 


7(, + 0.20 
■ '"-0.17 


-1 


= 9 + 0.16 
■"^-0.18 


0.484 




0.65355 


0.00427 


0, 


.95269 


32 


(13.20, 15.39) 




4. 


1 +0.7 
-^-0.6 


> 4.13 












1 


0.38152 


0.00574 


0, 


.57632 


17 


(13.15, 16.00) 




5. 


n+3.4 
"-2.2 


°-^^-2.14 


0, 


7Q + 0.32 
■ ' ^-0.25 


-1 


7C+0.28 
■ ' "-0.33 


0.234 




0.38152 


0.00427 


0, 


.57632 


15 


(13.20, 14.75) 




5. 


0+1.5 


> 4.81 












1+2 


0.40227 


0.00574 


0, 


.57632 


19 


(13.15, 16.00) 


io.5t|:? 


6. 


g+3.B 
°-2.3 




0, 


„n + 0.33 

.»S_p 27 


-1 


74 + 0.28 
■ '^-0.31 


0.049 




0.40227 


0.00427 


0, 


.57632 


17 


(13.20, 14.75) 




6. 


0+1.7 
'^-1.4 


> 5.64 












1 


0.78567 


0.65126 


0, 


.95269 


19 


(13.37, 16.00) 




3. 


g+l.S 

"-1.2 


6-35+^.;;^ 


0, 


(5„ + 0.27 
■"^-0.21 


-1 


oq + 0.24 
■■^■^-0.27 


0.268 




0.81814 


0.65126 


0, 


.95269 


12 


(13.40, 15.39) 


5.6;i:« 


2. 


7+0.8 
' -0.6 


> 2.48 












1+2 


0.81814 


0.65126 


0, 


.95269 


22 


(13.37, 16.00) 




4. 


7+2.1 
' -1.4 


- „o + 2.65 
' ■^■3-2.14 


0, 


7c;+0.28 
■ ' "^-0.23 


-1 


45 + 0.23 
■^"-0.25 


0.263 




0.81814 


0.65126 


0, 


.95269 


15 


(13.40, 15.39) 




3. 


1 +0.8 
-^-0.6 


> 3.10 













Samples by Instrument 



E140M 


0.06351 


0.00574 


0.08091 


7 


(13.15, 14.25) 


10.5±^^,« 


9 6+^"-' 

"■"-4.9 


5 65+-'^ '" 

"■""-2.91 


73+i-"S> 

"■ '"-0.51 


_2 i9+U-«« 


0.281 




0.06351 


0.00427 


0.08091 


7 


(13.20, 14.17) 


7 7+3.9 
' ■ ' -2.7 


7 Q + 3.6 
'■^-2.6 


> 3.74 








E230M 


0.73910 


0.24010 


0.95269 


34 


(13.14, 16.00) 




^■"-1.1 




70+0.22 
"■ '"-0.18 


, 00 + 0. 17 


0.160 




0.73910 


0.24010 


0.95269 


25 


(13.25, 15.39) 


""-1.2 


3 fi + 0-7 
■^■"-0.6 


> 4.32 









Note. — Parameters from the maximum-likelihood analysis for /(A'"{C"^^)) — fci4 (A'"(C"'"^)/No)"-'^ , where log No — 14. The CIV doublets were 
divided into several sub-samples by the group G and the redshift range. For each sub-sample, the first line refers to the maximum-likelihood analysis 
and the second line, to the observed quantities. dA/civ/dX, listed in the first sub-sample row, is the integral of /(7V(C + ^)) from logA^(C+^) - 13 
to infinity with the best-fit ki4 and a. jv (see Equation |18| l . Also in the first sub-sample row, the integrated dA'c iv /^^ = dAc iv /dX ■ dX/dz, where 
the latter term is the derivative of Equation[9] evaluated at {z). The observed dA/civ/dz and dA'civ/dX are from the sum of the total number of 
doublets, weighted by the pathlcngth available to detect the doublet (based on its Ar(C'^^); see Equation |17| l. 0^+3, listed in the first sub-sample 
row, is the integral of /(A^(C+^)) ■ 7V(C'^^) from 13 < log7V(C+^) < 15 with the best-fit /ci4 and an (sec Equation[22j. The observed ft(~.+3 were 
from the sum of the unsaturated doublets, as given by J\f (sec Equation |20| l . Pks is the significance of the onc-sidcd K-S statistic of the best-fit 
power law. 



Table 7 

/(W^r,i548) AND (Wciv/clX SUMMARY 



(1) 




(2) 


(3) 




(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


G 






2; 








(mA) 


dA^civ/dz 


dA/'civ/dX 


fe3 

(mA-1) 


aw 


Pks 


















Samples 


by Rodshift 










1 


0, 


.65355 


0.00427 


0, 


.95269 


38 


(47, 2423) 


yg+i.'J 

' ■'J-i.e 


, 7+1. U 
■5- ' -0.8 


-^■O''-0.40 


l.DU_Q 26 


0.714 
















(53, 2363) 


5 7+11 

°- ' -0.9 


3 4+° '' 




1+2 


0, 


.65355 


0.00427 


0, 


.95269 


43 


(37, 2423) 


7.8+=,-? 


4 1+1- = 


1 7O + 0.53 
-'-■'8-0.42 


, t-7 + 0.22 
-^■O'-0.23 


0.374 
















(52, 2363) 


6.7ti;? 


o Q+0.7 
■^■^-0.6 




1 


0, 


.38152 


0.00427 


0, 


.57632 


18 


(47, 896) 
(53, 852) 


7 c + 11.5 
'■0-2.4 

6 2+^ '' 


4 8+''-^ 
4 6+^* 


-0.88 


-1 60+°''** 
^■"'^-0.49 


0.200 


1+2 


0, 


.40227 


0.00427 


0, 


.57632 


20 


(37, 896) 


8 4+"-° 

°-*-2.8 


5.3+';? 


■^■^"-0.91 


-1 58+°'"' 
J--OO_0.45 


0.038 
















(52, 852) 


7 7+2.3 
' ■ ' -1.9 


5 7+1-7 
0- ' -1.4 




1 


0, 


.81814 


0.65126 


0, 


.95269 


20 


(61, 2423) 
(67, 2363) 


7 6+^" 

5 5+^-^ 
0-0-1.2 


3.6+t? 

o 7+0.8 
^■'-0.6 




-1 46+°-^' 


0.501 


1+2 


0, 


.81814 


0.65126 


0, 


.95269 


23 


(61, 2423) 
(67, 2363) 


6.3llj 


, 1+0.8 
■^■-^-0.6 




_1 4g+0 32 

^■^■^-0.34 


0.624 



Samples by Instrument 



E140M 


0.06351 


0.00427 


0.08091 


8 


(47, 315) 
(53, 264) 


Si 7.4 


fa 6.7 

6 1+^* 

"■^-2.3 


2 36+"-^-* 


-1 89+'''-* 


0.248 


E230M 


0.73910 


0.24010 


0.95269 


35 


(37, 2423) 
(52, 852) 


5.211:? 


4 3+2 = 

9 7+0.7 
' -0.6 


1 qo+°-=^ 

.^■'"^-0.48 


-1 42+°-24 


0.355 



Note. — Parameters from the maximum-likelihood analysis for f{Wr) — (Wr /Wr ,0)°^^ , where Wr .0 — 400 mA. The tabulated 
information is similar to that presented in Table [ol except that the integrated line density limit is Wr — 50 mA. 



Table 8 

/(iV(C+3)) AND /(Wr,i548) CORRELATION MATRICES BY ABSORBER 









/(Ar(C+^)) 






/(W,,1548) 




G 


z 


(fel4) 


(on) 




(fc3) 


(aw) 


rki,ai 






(cm2) 






(mA-') 






1 


< 1 


0.65 ± 0.04 


-1.50 ± 0.12 


-0.63 


1.57 + 0.62 


-1.59 ± 0.33 


0.97 


1+2 


< 1 


0.75 ± 0.04 


-1.52 ± 0.12 


-0.69 


1.76 ± 0.72 


-1.56 ± 0.33 


0.97 


1 


< 0.6 


0.75 ± 0.05 


-1.76 ± 0.21 


0.43 


1.96 ± 0.49 


-1.60 ± 0.34 


0.99 


1+2 


< 0.6 


0.84 ± 0.05 


-1.74 ± 0.22 


0.33 


2.18 ± 0.51 


-1.58 ± 0.31 


0.99 


1 


> 0.6 


0.60 ± 0.11 


-1.40 ± 0.17 


-0.78 


1.57 + 0.66 


-1.44 ± 0.38 


0.87 


1+2 


> 0.6 


0.72 ± 0.13 


-1.46 ± 0.16 


-0.76 


1.79 ± 0.81 


-1.47 ± 0.41 


0.88 



Note. — Mean values and correlation coefficient from jack-knife analysis of power-law fit to 
f (N {C+^)) and f{Wr). Each absorber was iteratively excluded from the maximum likelihood analysis 
of the frequency distributions. 
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Table 9 

/(iV(C+3)) AND /(Wr,i548) CORRELATION MATRICES BY SIGHTLINE 















/(W,,154S) 




G 


z 




("«) 






(aw) 








(cm^) 












1 


< 1 


0.63 + 0.11 


-1.50 ± 0.10 


-0.12 


1.54 + 0.62 


-1.58 + 0.31 


0.88 


1+2 


< 1 


0.73 ± 0.11 


-1.52 ± 0.10 


-0.14 


1.73 + 0.67 


-1.55 + 0.32 


0.92 


1 


< 0.6 


0.72 ± 0.20 


-1.76 ± 0.27 


0.70 


1.89 + 0.80 


-1.61 + 0.39 


0.94 


1+2 


< 0.6 


0.81 ± 0.20 


-1.75 ± 0.27 


0.45 


2.11 + 0.79 


-1.58 + 0.37 


0.91 


1 


> 0.6 


0.56 ± 0.11 


-1.40 ± 0.18 


-0.72 


1.52 + 0.60 


-1.43 + 0.37 


0.82 


1+2 


> 0.6 


0.69 ± 0.09 


-1.46 ± 0.14 


-0.39 


1.75 + 0.65 


-1.46 ± 0.41 


0.90 



Note. — Mean values and correlation coefficient from jack-knife analysis of power-law fit to 
/{A^(C+'^)) and fiWr). All absorb crs from each sightlinc were iterativcly excluded from the maxi- 
mum Ukehhood analysis of the frequency distributions. 



where iVo = lO^'^cm"^ and Wrfi = 400mA and 
the subscripts on the coefficients k indicate the 
normahzation {e.g., ki4 — fc/IO^^cm^). We used 
the conjugate gradient method to maximize the 
hfcchhood function £ and to simultaneously fit for 
the coefficient k and exponent a. From the 1-a 



error ellipse where In £ — In C„ 



> - 



1.150 



we 



estimated the errors in k and a (see Figure [9]). 

We de rived £ in a similar manner to that out- 
lined in Storrie-Lombardi et al. |l99d) and de- 
tailed in our Appendix |B] The main difference 
was to include the observed number of strong, 
saturated absorbers, where we only had a lower- 
limit estimate of A^(C+'^), as a constraint in £. 
The best-fit power law must allow for a number 
of strong (saturated) absorbers consistent with 
the observed number. We set the saturation 
limit to log iV(C"'"'^) = 14.3, which we determined 
empirically (see Figure |4]). In the maximum- 
likelihood analysis, the unsaturated doublets with 
logA^(C"'"^) > 14.3 were counted as saturated ab- 
sorber, and the two saturated (G = 1-1-2) doublets 
with logAf(C+^) < 14.3 were excluded. Without 
the saturation term in C, the best-fit a would have 
been significantly steeper, depending on the inte- 
gration limits, since it would have been a strong 
statistical statement to not detect any strong ab- 
sorbers. 

For the same reason, the choice of integration 



^^We empirically determined the In £ — ln£max > —1.15 con- 
straint. This limit defines a contour that contains 68.3% 
of the likelihood surface. For a Gaussian distribution, the 
l-cr contour would be defined by In £ — In^max > —0.5 



limits in the maximum-likelihood analysis influ- 
enced the result. We set the lower integration limit 
to the smallest observed value less one sigma for 
the sample analyzed, e.g., Nmin — cr^,min- For 
/(7V(C+^)), the saturation limit was 10 cm 
and the upper limit "infinity" was lO^^cm"^. For 
/(W^r,i548)j the upper limit was the largest ob- 
served value plus one sigma Wr,max + cfWr-max- 

The best-fit power-law f{N{C+^)) for the G = 
1 and the G = 1-1-2 samples are shown in FigurefTUl 
and enumerated in Table [6l For the results of the 
maximum likelihood analysis of /(Wr,i548)j see 
Table [T] They were consistent within the errors. 
To examine the temporal evolution of C IV ab- 
sorbers in our survey, we divided the G = 1 sample 
into two redshift bins, defined by, approximately, 
the median redshift: Z1548 < 0.6 and 0.6 < z < 1 
(see Figure [TT|) . 

We performed a jackknife re-sampling analysis 
to estimate the errors of and correlation between k 
and a from our maximum-likelihood analyses. For 
the jackknife, each C IV doublet i was excluded 
from the full sample of A/" absorbers. Then, the fre- 
quency distributions /(A(C+^)) and /(Wr, 1548) 
were re-fit with integration limits set as described 
previously. This was done for the G = 1 and G — 
1+2 samples for all three redshift cuts (^1548 < 1, 
< 0.6, and > 0.6; see Table[S]). The variance in the 
e.g., coefficient distribution ki from the jackknife 
is: 



AA- 



{k))\ 



(15) 
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where (k) is the mean of the best-fit values from 
the re-sampling analysis. The same applies to the 
exponent distribution ai. The mean coefficients 
and exponents from the maximum-likelihood anal- 
yses to the various redshift samples agreed with 
the results from the analyses of the full samples 
{i.e., compare Tables [6] and [7] with Table [8|). 

The correlation between the parameters ki and 
at is the ratio of the covariance to the product of 
and CTq, . : 



(16) 



The coefficient and exponent were strongly cor- 
related for the /(W^r,i548) power-law fits (i.e., 
Tki.oi ~ 1). By accounting for the saturated dou- 
blets in the maximum-likelihood function, kn and 
aN were less correlated for the /(A^(C+^)) fits. 



Sin ce C IV absorbers cluster ()Fernandez-Soto et al 
19961) . we also performed a jackknife of the sight- 



line and iteratively fit the frequency distributions 
after excluding all absorbers from each sightline. 
The results are compiled in Table [HI they agree 
well with results for the full sample (Tables [6] and 
[7]). We did not re-estimate AX without the ex- 
cluded sightline since each sightline contributes 
< 10% to AX of any redshift range. 

5.2. C IV Absorber Line Density 

We measured the density of C IV doublets with 
a minimum column density as follows: 

djVciv 



dX 



Ni>Nii, 



1 

AX{N, 



(17) 



'dJ\f/dX 



E ( 



Ni>N,i 



AXiNif 



The absorber line density can be measured for any 
equivalent width limit in a similar fashion. We es- 
timated the contribution of the low-number statis- 
tics in Af by adding the Poisson counting variance 



to cr, 



In Figure fT2l we show cLA/civ/dX for 
the full G = 1 and G — 1+2 samples as functions of 
column density (top panel) and equivalent width 
(bottom panel). The weakest system we found 
had logiV(C+3) = 13.2 and W^,i548 = 52 mA. At 
this limit, dAbiv/dX = i.ltol for the G = 1-^2 
sample. 
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Fig. 9.— /(7V(C+3)) and /(Wr,i548) 1-<J ellipses 
from the maximum- likelihood analysis of the G 
— 1 sample. The sample was divided by redshift: 
z < I (black ellipse and points), z < 0.6 (light 
gray, dashed), and 0.6 < z < 1 (gray, dotted). 
The plus signs indicate the best-fit values (see Ta- 
bles m and [T]) . The diamonds indicate the coef- 
ficient and exponent that define the 1-a error on 
the integrated and/or dA/civ/dX. The co- 

efficient kii and the exponent aN were not highly 
correlated for /(A^(C"'"'^)) (top panel) due to the 
treatment of saturated absorbers in the maximum- 
likelihood analysis (see Appendix |B]). The val- 
ues were highly correlated for f{Wr,i54s) (bottom 
panel). 

In Figure [121 we extrapolated our dAciv/dA 
to lower A^(C"'"'^) and Wr by integrating the power- 
law fits of the frequency distribution over e.g., 
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Fig. 10. — Column density frequency distribution 
for the full z < 1 sample. The G = 1 observa- 
tions are the black diamonds, and the long-dash 
line indicates the best power-law fit. The observed 
and fitted G = 1+2 f{N{C+^)) agrees very weh. 
The AODM column densities are lower limits for 
saturated systems; our adopted saturation limit 
logA^(C+^) = 14.3 is indicated (vertical, dotted 
line) . 
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(18) 

The summed and integrated C IV absorber line 
density by A^(G+'^) and Wr,i548 are listed in Ta- 
bles [H] and [71 respectively. The error in dAc iv / 
from the power-law model of /(A^(G+^)) was 
defined by the error ellipse in the maximum- 
likelihood analysis. However, the extrema in 
dAfcw /dX allowed by the 1-cr eUipse did not 
occur at the extrema of fci4 and aN- As shown 
in Figure [9] (bottom panel), the extrema in the 
1-CT-allowed values of dAfcw /dX occurred on the 
sides of the ellipse where fci4 and apf were max- 
imized (or minimized). Again, the equations are 
similar for line density as a function of equivalent 
width. 

5.3. Mass Density 

In principle, ri^+a is the ratio of the mass in 
ions relative to the critical density pc,o- In 
practice, the observations limit Cl^+a to include 



Fig. 11. — Column density frequency distribution 
for two redshift bins of the G = 1 sample. The 
best-fit UN for the 21543 < 0.6 f{N{C+^)) (black, 
long-dash line) is steeper than that for the 0.6 < 
•2^1548 < 1 bin (gray, dashed line). As in Figure 
I10| the best-fit aj^ and kn from the maximum- 
likelihood analysis are shown, as is the saturation 
limit log7V(C+^) = 14.3 (vertical, dotted line). 



only C IV absorbers within a range of column 
densities. The lower A^(C"'"'^) bound reflects the 
limit where the observations can confidently de- 
tect and identify C IV doublets, typically Nmin ~ 
lO^^cm"^. The upper bound is usually Nmax = 
10^5 cm-2. Doublets with log7V(C+3) > 15 are 
rare and often associated with galaxies (for exam- 
ple, the DLA C IV doublet towards PG1206-)-459 
with log7V(C+^) > 15.4). In addition, there has 
been no observed break in /(A^(C+'^)), so the in- 
tegrated is infinite. 

The metallicity and ionizing background of the 
intergalactic medium affect the C"*"^ mass den- 
sity, and the observed evolution of over time 
(or redshift) add a constraint to the changes in 
the cosmic metallicity and ionizing background. 
While one expects that the former increases mono- 
tonically with time, the latter is believed to be 
decreasing since z k, 1. VLq+'s is the integrated 
column density- "weighted" frequency distribution 
/(7V(C+3)): 
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Fig. 12. — Absorber line density as function of 
iV(C+3) and Wr,i548- The summed dAfcw/dX 
from the G = 1 and 1-1-2 groups are shown with 
the sohd black and gray crosses, respectively (see 
Equation [T7)) . The dashed lines show the inte- 
grated cWc w/dX from the fits to f(N(C+^)) and 
/(W^r,i548) (see Equation [TH]). 



where Hq = 70kms~^ Mpc^^ is the Hubble con- 
stant today; mc = 2 x 10~^^g is the mass of 
the carbon atom; c is the speed of light; and 
p^ Q = 3i?^(87rG)-i = 9.26 x IQ-^Ogcm-^ for our 
assumed Hubble constant. 

The observed f2c+3 can be approximated by the 
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Fig. 13. — C^"^ mass density relative to the 
critical density over age of Universe. The in- 
tegrated (see Equation for the G = 
1 sample, divided by r edshift, are the (orange) 



stars. The value from lOanforth fc Shulll (|2008h 



is the (green) triangle and is not an indepen- 
dent measurement of ilc+a . The grey boxes in- 
dicate the binned, cosmology-a djusted ^n + s from 
several high -redshift 
[Pctt ini et al 



studies: iSongaila 



Sc aiinapieco et al]~~(l2006l) 



( 20031): iBoksenberg et al 



and iRvan- Weber et al 



Becker et al 




" pooi l (see Appendix p. 
The horizontal extent of the boxes indicate the bin 
size. The vertical extent represents the maximum 
range spanned by flQ+i ± aa (the published er- 
rors were adjusted to be la, as necessary). The 
average time and error-weighted Qc+3 per bin are 
shown with black asterisks. 



sum of the detected C IV absorbers: 
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(jLanzetta et al.lll991r ). We list the summed il^+a 
values for various redshift samples in Table [H We 
only include absorbers that are unsaturated in at 
least one line, typically with log A^(C-'"'^) < 14.3, 
because we only have a lower limit to iV(C-^'^) for 
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saturated absorbers. Therefore, the summed f2c+3 
is a lower hmit, since some saturated absorbers 
may have logiV(C+^) < 15 and could have been 
included in the sum. 

To measure a value for ri^+a , we use our best-fit 
/(iV(C+^)) power law fEquation[T3)) and integrate 
Equation [19] analytically: 



C+3 



CPcfl 



UN 



N. 



2+qn 



(22) 

With this, we can measure fl^+a over the column 
density range of unsaturated and saturated ab- 
sorbers or even extrapolate the model to different 
ranges. Like the errors for cLA/civ/dX (see § 15.21 
the errors for the integrated flc+s were taken as 
the extrema of the 1-cr /(A^(C"'"'^)) error ellipses 
(see top panel. Figure [H). 

The choice for the limits of integration Nmin 
to Nmax was critical, since there has been no ob- 
served downturn in /(iV(C^'^)), measured at any 
redshift. We used 13 < \ogN{C+^) < 15 in our 
analysis (see Table |6] and Figures [13] and [TS]) be- 
cause this range reflects that of the observed col- 
umn densities and overlapped with the majority 
of st udies used for comparison {e.g., S ongaila 
l200l[ ). The integrated Q^+a is most influenced 
by the strongest absorbers, because they con- 
tain substantially more mass. For example, ac- 
cording to Equation [22] and with ajv = —1.5, 
the 14.3 < logAf(C+^) < 15 component con- 
tributes « 60% of flc+s , integrated over the range 
13 < log A^(C"'"'^) < 15. In comparison, the 
13 < logA''(C+'^) < 13.7 component accounts for 
« 14%. The proportions change as the power- 
law exponent changes, increasing the lower col- 
umn density contribution as aN decreases. A de- 
tailed discussion of the effects of the column den- 
sity range on the summed flc+s is given in Ap- 
pendix [C] 

6. Discussion 



6.1. Comparisons with Previous Results 

In lDanforth fc Shulll ^00^), (Wcw/dz 
for Wr > 30 mA. They measured the absorber 
line density by summing their observed number 
of doublets (24) and dividing by the unblocked 
redshift pathlength {Az = 2.42). They esti- 
mated Az by identifying regions of the spectra 



where the C IV 1548 line could be detected at 
Wr > ^cTWr 1550 at > 2(7^^- We measure 
dNciY/Az = 6.21^;^ for the z < 6.6, G = 1 dou- 
blets with 1^^,1548 > 53 mA (see Table [l]). The 
> la difference between our dAfciy/dz values is 
largely due to the equivalent width limit. For 
Wr > 50 mA, dA/biv/dz = itl (C. Danforth, pri- 
vate communication), which is within 0.3cr of our 
value. Another source for the discrepancy is the 
detection criterion used {e.g., both lines detected 
with Wr > 3(7 ) , which would affect the number 
of doublets included and the unblocked redshift 
pathlengthF^ 

The C IV absorber line density has not changed 
significan tly since z = 5 for log iV(C^^) > 13 ab- 
sorber s. Songaila ( 200ll ) and iBoksenberg et al. 
(l2003i) did not detect any redshift evolution for 
log A^(C+3) > 13 doublets over 1.5 < z < 4.5; they 
measured dJVcw/'iX w 3. In Figure [T^ we com- 
pare t he cLA/civ/dA^ measurements from Songailal 
(|200l[ ) with the current study. In order to com- 
pare <Wcw/dX for logiV(C+3) > 13.2, we use 
the best-fit coefficient (fci4 = .63 cm^) a nd ex po- 
nent {on = -1.8 ± 0.1) fromlSongaiy (|200ll ) in 
Equation[18] <Wcw/dX = 1.9±0.2 (adjusted for 
cosmology). For the current survey, the summed 
absorber hue density is dA/civ/dAT — S.A'^q^ at 
(z) = 0.654. This is a 1.8 ± 0.4 i ncrease over th e 
integrated dAfcw/dX value from' Songailal (2001). 
However, since she did not provide an error for 
the coefficient, we likely underestimate her error 
on dA/civ/dX above. 

The best-fit = -l-^^toHi for the lower- 
redshift bin was consistent with a ^ = —1.79 ± 
0.17 from lDanforth fc Shulll (|2008h . which covered 
Z1548 < 0.12. 

There has been no consensus on the be st power- 



law ex ponen t for /(jV(C+^)) a t z > 1. ISongaila 
(|200ll ) and iBoksenberg et all (|2003h conducted 
surveys similar to the current study, and they mea- 
sured aN = -1.8 ± 0.1 (2.9 < z < 3.54) and 
aN = —1.6 (1.6 < z < 4.4), respectively. For 
our 0.6 < Z1548 < 1 sample, t he best-fit a^ = 



1.39 ln '^7 agreed better with iBoksenberg et al 



(|2003[ ) but was consistent with both studies at the 
2(7 level. 

More detailed comparison with other studies of 
C IV absorbers at z < 1 are given in Appendix |D] 



i^We measure Az = 2.2 for the E140M spectra. 
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Fig. 14. — cWciv/d-'^ as a function of red- 
shift. The number density of C IV absorbers 
has not changed significantly since z — 5. The 
(orange) stars are our summed z < 1 values 
for logiV(C+3) > 13.2 (see Equation [HI) . The 
solid and dashed horizontal lines are the inte- 
grated dAfcw /dX and l-t r error, r espec tively, 
from the best-fit f{N) from ISongaili (|200l[) . She 
only fit her 2.9 < z < 3.54 observations, but 
the result agreed with the full 1.5 < z < 5 
sample. For reference, dJ\fmv / dX = Af/AX 



from the tabulated M and A X inlSongaila 



12 < logiV(C+3) < 14.9) and lPettini et al 



2001 



2003 . 



12.5 < logiV(C+^) < 14) are shown with the gray 
pluses and cross, respectively; the error bars as- 
sume Poisson counting statistics on Af. We have 
adjusted AX for differences in cosmology (see Ap- 
pendix [C| . 



6.2. nr 



Evolution 



We have measured a statistically significant in- 
crease in il(j+3 at z < 1 compared to the roughly 
constant value observed at 1 < z < 5. The error- 
weighted average of the 1 < z < 5 measurements 
in Figure [13] is Hc+s = (2.2 ± 0.2) x 10"^ at z 
3.240. Our integrated Oc+3 = (6.20^^;^) x 10"^ 
at (z) = 0.654 is a 2.8 ± 0.7 increase. We rec- 
ognize the significance of the increase in Qq+3 
at z < 1 because we have carefully adjusted the 
1 < z < 5 values of other aut hors (see Appendix 
Danforth fc Shulll (120081) stated that their 



ED 

Qr,+3 value agreed with that of lScannapieco et al.l 
( 20061 ) without considering the effect of the col- 



umn density limits. iFrve et al.l (|2003l ) interpreted 
their z < 0.1 flc+^ value as an increase but with- 
out comment on its significance. 

To estimate the rate of evolution of the 
mass density, we adopt a simple linear model: 
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flo + aitan 



(23) 
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age ai 



2ta,eCOVoi, 



where tage is the age of the Universe and COVqi 
is the covariance of the intercept ao and the 
slope ai. The observed 1-a uncertainties in 
rip+awere used in the x^-minimization algo- 
rithm. The results of the linear regression for 
the 1 < z < 5 (IGyr < tage < 6Gyr) and z < 5 
{tage > 1 Gyr) samplcs are shown in Figure 1151 
We included the following Qq+s measurements, 
adjusted for cosmology and /or co l umn density 
range (see A ppendix ICl): Songailal ( 200lL 1.5 < 
z < 4.5) ^1 IPettini et al.l (120031. as revision of 
Songailal (12001) z ~ 4.7 value); Boksenberg et al 



( 20031): IScannapieco et al.l (|2006[ ): and the current 



studyo 

Whether or how flQ+a evolves from z = 5 ^ 1 is 
not statistically constrained, assuming the simple 
linear model. The best-fit slope for this range was 
ai = (0.15 ±0.3) X lO^^Gyr"^ with the intercept 
ao = (1.9 ± 0.6) X 10-« (at tage = OGyr) and 
probability P^2 = 34%. 

The linear regression for the z < 5 data in- 
dicates a statistically significant trend in Q^+a 
evolution. The fitted parameters are consistent 
with evolving shghtly: ai = (0.42 ± 0.2) x 

10-8 Gyr"^ and intercept ao = (1.33 ±0.5) x 10"* 
(Pv2 = 36%). Several high-redshift studies 



( Songailal200lt|Pettini et al.l2003l: [Boksenberg et al 
20031 ) concluded that ^c+a evolved very little to 



not at all for z « 5 — *■ 1, which is consistent with 
the previous linear regression of the 1 < z < 5 
data. Incorporating the new z < I measurements, 
there is evidence that f^c+s has slowly but steadily 
increased since z ~ 5, at the 97% confidence level. 

We acknowledge that our accounting for differ- 
ences in cosmology and column density range were 



I'^The results from IPettini et al] ||2003| ) are considered a re- 
vision to the 4.5 < 2 < 5 bin since [Pettini ct al. ( 2 00^ 
used higher S/N data. The 5 < z < 5.5 bin from Songailal 
(2001), which only included one detected absorber, was ex- 
cluded. 

l-*We excluded the iDanforth fc Shulll l!200^ measurement 
since our z < 0.6 measurement includes the E140M data. 
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imperfect. However, whether or not we adjust for 
differences in the iV(C+'^) range does not signifi- 
cantly affect the previous results: our fJ^+a value 
is still a statistically significant increase over fig+a , 
and the rate of evolution (oi) from z = 5 — > is 
also significant, assuming the simple linear model. 

We emphasize that Equation [221 was not phys- 
ically motivated. The temporal evolution of 
f2c+3 is influenced by multiple, complex phys- 
ical processes {e.g., star formation, UV back- 
ground). In addition, only C IV absorbers with 
13 < logAf(C+'^) < 15 are included in the mea- 
surements of 5^c+3. However, as we discuss in § 
16.31 the physical nature of doublets in this column 
density range likely changes over the 12 Gyr from 
z = 5 — > 0. The linear regression analysis was 
simply a secondary way to gauge whether our ob- 
served flQ+3 values suggest a significant increase 
compared to the z > 1 values. In the context of 
this simple model, we rule out null evolution at 
high confidence. 

6.3. Changing Nature of C IV Absorbers 

In § 16.11 we noted that dA/civ/dX for 
log A''(C"'"'^) > 13 has not evolved with redshift 
(see Figure [Til) . However, there is evidence for 
evolution of dAfcw /dX when only the strong sys- 
tems (Wr,i548 > 150 mA an d loKiV(G+3) > 13.7 



1 2 3 5 



are included (ISteidell il990l : iMisawa et alJIioQ- 



i 



Boksenberg et all 2003 ) . The weak systems dom- 
inate the number counts, resulting in no evolu- 
tion of (Wcw/dX for the full sample. Since the 
stronger doublets dominate the C"*"^ mass den- 
sity, we observe the increased number of strong 
C IV doublets in the significant increase of 17^+3 
at z < 1. 

Even without limiting redshi ft-evolution analy- 



Boksenberg et al 



sis to the strongest absorbers 

(|2003') noted that complex C IV systems have 
higher mean column densities at lower redshift, 
where "complex" indicates that t he number of 
Voigt profile components is > 7. iMisawa et al 



(|2002[ ) predicted this result by tying their ob- 
served inc rease of high-W^ doub l ets at low red- 
shift with iPetitiean fc Bergeronl (| 19941 ) observa- 
tions that higher Wr systems have more compo- 
nents. From visual inspection of the W, 



r,1548 



> 



150 mA doublets, we see that more than half have 
multiple, prominent components. 
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Fig. 15. — Linear regressions of flc+s over age 
of Universe. The (red) dashed lines are the x^- 
minimization fit to subsets of the data described 
in Figure [131 and the (red) dotted lines indicate 
the l-cr range of the model. In the upper panel, 
the slope for the 1 < z < 5 fit is consistent with 
no evolution ai = (0.15 ± 0.3) x 10"^Gyr~^ For 
z < 5, in the lower panel, the slope is ai — (0.42 ± 
0.2) X lO^'^Gyr"^ The evidence that the 0+^ 
mass density has been slowly increasing since « 
1 Gyr requires the measurements at low redshift. 



Perhaps the trend to higher column den- 
sity, more multi-component C IV absorbers at 
low redshift supports the claim that ~ 50% of 
logA^(C+'^) > 13. 3 doublets are associated with 
galactic outflows ( Songailal 20061 ). Alternatively, 
the strong absorbers might trace infall or high- 
velocity cloud-like halo gas. In either case, the 
complex profiles might be a result of the het- 
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erogeneous nature of outflowing or infalling gas. 
Observations at low-redshift have indicated that 
C IV absorption often reside in ga laxy halos, on 
« 100 kpc scales (|Chen et al.ll200l[ ). 

What C IV absorption traces at all redshifts af- 
fects the interpretation of the evolution of fl^+a . 
In § 16.21 we discussed the time evolution of ^Iq+s 
for 13 < logiV(C+3) < 15 doublets at z < 5, and 
we showed that there has been a significant in- 
crease at z < 1 and that null evolution is ruled out. 
However, if the majority of 13 < logiV(C+^) < 15 
C IV doublets are (low-density) intergalactic at 
z > 1 while the majority are circum-galactic at 
z < 1, then the Qc+^ observations at high and low 
redshift are not directly comparable. In a future 
paper, we will explore the changing nature of C IV 
doublets by comparing the properties and environ- 
ments of C IV absomtion in cosmological hydro- 
dynamic simu lationa^^l with observations (also see 
Cooksevl[2009h . 



geted bright quasars, without prior knowledge of 
e . (?. , Mg II, damped Lyg , or Lyman limit s ystems 



6.4. Sightline Selection Bias 

Several sightlines in our sample were originally 
observed for reasons that increase the likelihood 
that C IV absorption would be detected. For 
example, there were quasars targeted because in- 
tergalactic Mg II absorption was observed in op- 
tical spectra. Studies have shown that Mg II 
absorbers frequently exhibit C IV absorption 



(jChurchill et al.lll999a[ ). Including such sightlines 



in our survey might have increased our detected 
doublets and biased our results to have more C IV 
doublets than would be observed in an wnbiased 
survey of the IGM. On the other hand, the sen- 
sitivity limits of high-resolution spectrometers on 
HST (prior to the installation of the Cosmic Ori- 
gins Spectrograph) imply only a small number of 
quasars could have been observed, independent of 
known absorbers. 

To explore the bias in our survey due to 
the sightline selection criteria, we compared 
our C IV absorber line density cWciv/dz with 
Barlow Sz Tvtler ( IQQSh . who performed a sur- 
vey for strong C IV absorbers in 15 sight- 
lines from the HST Faint Object Spectrograph 
(FOS) archives at the redshifts of known Lya 
absorbers. The FOS quasar key line project tar- 



(|Bahcall et all Il993h . iBarlow fc Tvtleil (|l998h 
measured dA/biv/dz = 2.3 ±0.9 for Wr > 400 mA 
and 0.2 < 2:1548 < 0.8 (six doublets). 

We d etected six G — 1 C IV doublets that 
met the Barlow fc Tvtleil ( 19981 ) criteria and find 
for this sub-sample, dAciv/dz = 1.3to;5. Thus, 
our measured absorber line density agreed with 
the unbiased result within la. Despite including 
sightlines observed with prior knowledge about the 
IGM, we agreed with an unbiased survey and, in 
fact, record a smaller incidence of strong C IV ab- 
sorbers. 

As a second check for sightline bias, we reana- 
lyzed our survey after excluding 11 G = 1 C IV 
doublets that were possibly associated with tar- 
geted Mg II absorbers, DLAs, or LLS (see Table 
10l)I!!| For this sub-sample, dA/biv/dz = 0.4^°;^ 
for the two G = 1 C IV doublets with Wr > 
400mA and 0.2 < z-]mh < . 8, wh ich under- 
estimates the Barlow fc Tvtler ( 19981 ) results by 
1.8(7. The best- fit coefficient and exponent from 
the maximum-likelihood analysis of /(Wr,i548) 



were k-t 



1.40 



+0.67 , 
-0.51 



lA and aw — —1.51 



+0.39 

-0.41' 



which lie within 0.5tT of the values in Table [T] 
For /(7V(C+^)), the best-fit values were fci4 = 
0A9toltcm^ and ajv = -l-7lto;26, which were 
within 1(7 of the values in Table [H respectively. 
The change in coefficient and exponent change the 
integrated ^c+s = (SMti lt) ^ 10"^' lower than 
that from the full sample but a 1.8 ± 0.5 increase 
over the error-weighted average of the 1 < z < 5 
measurements. 

Excluding the 11 C IV doublets possibly asso- 
ciated with known absorption features yields an 
integrated Qc+^ within 3(7 of the value extrap- 
olated from the model to 1 < z < 5 (see top 
panel, Figure I15p . However, by excluding the 
11 absorbers, most with log7V(C''"'^) > 14, we 
under-sampled the strong doublets as expected 



^^The simulations are from the Overw helmingly Large Sim- 
ulations project dSchave et al.|[2"009D . 



"Several sightlines were targete(i for a specific absorp- 
tion system but the archival spectra did not cover the 
C IV doublet. The following are the sightlines and 
bias: PKS0454-22. Mg I I at Zabs = 0.6248 and 0.9315 
dChurchill fc Le Brunll998l '): HE0515-4414, DLA at z^bs = 
1.15; MARK132, LLS at z^,,, = 1.7306; and PKS1127- 
145, DLA at z^bs = 0.312. The PG1634+706 sightUne 
was also targeted for a D/H study and Mg II absorber at 
^abs = 0.9902, but we did not include that echelle order in 
our co-added spectra. 
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from the iBarlow fc Tvtlej (|l998[) study, and the 
high- A'' (C"'"'^) systems dominate the C"'"'^ mass den- 
sity. In the attempt to eliminate the effects of the 
sighthne selection bias, we introduced a possible 
bias against strong C IV doublets. We have pro- 
ceeded, therefore, with the full sample of archival 
spectra. 

6.5. C IV Doublets without Lya Absorp- 
tion 



U 

^ 14.0 
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Fig. 16. — C"*"^ and H I column densities. The 
symbols indicate the location or limit of all C IV 
doublets with both lines detected at Wr > iuw^. 
and with Lya coverage. When possible, the H I 
column densities included the constraint given by 
Ly/3 detection. The G = 1 sample are highlighted 
by the filled black square. The upper limits are 
2ct. 

We conducted a blind survey for C IV doublets, 
relying initially on the characteristic wavelength 
separation of the doublet. From the resulting list 
of candidates, we drew on other known charac- 
teristics of C IV doublets to distinguish the final 
sample (see §|3]and Table [3]). Although the detec- 
tion of associated Lya absorption was a dia gnostic 



and based on previous ob servations (jEUison et al 



19991 : ISimcoe et al.l it was not required. 



C IV doublets withou t Lya absorption have been 
observed at z > 1 ( Schave et al.l 120071 ). Here we 
discuss the likelihood that some C IV doublets in 
our sample were actually due to Lya forest con- 
tamination (see § 14. 3p and if we indeed detect any 
"naked" C IV absorption. 



Our sample includes two G = 2 C IV dou- 
blets without Lya absorption detected (A^hi < 
lO"'^^ '' cm~^) that were in the Lya forest: 21548 = 
0.83662, logiV(C+3) = 14.2 doublet toward 
HS0747-H4259 and Z1548 = 0.87687, logiV(C+3) « 
14 doublet in toward HS0810-h2554. In the scat- 
ter plot of A(C+3) versus A^hi, the HS0747-h4259 
system has the highest C"'"'^ column density for 
any logA'^ni < 14 system (see Figure [16]). The 
HS0810-I-2554 system does not have a useful con- 
straint on A'^Hi (see Table [5]) and could not be 
plotted. Though these two systems c ould be high- 
metallicity systems akin to those in ISchave et al 



( 20071 ) , it is more likely that their atypical column 
densities are due to their low signal-to- noise ratios, 
which increases the uncertainties in continuum fit- 
ting and column densities. The HS0747-I-4259 and 
HS0810+2554 E230M spectra have S/N = 6pix~^ 
and 3pix^^, respectively (see Table [T]). 

The low S/N also increases the likelihood that 
Lya forest lines will be (incorrectly) identified as 
C IV doublets (see § 14. 3[) . In the most extreme 
scenario, the maximum absorption pathlength in 
the Lya forest is AA = 8.9 for the strongest 
Monte-Carlo absorbers (logA(C+^) > 14.1 or 
W'r,i548 > 157mA). The estimated rate of Lya 
forest lines masquerading as C IV doublets was 
dAA/dA < 0.11 for 0.8 < z < 1.5 (see Figure 
[8|). Therefore, the maximum expected number of 
C IV- like Lya pairs would be 1^^, which is con- 
sistent with Lya forest lines coinciding to be the 
two G = 2 doublets in question. 

Next, we evaluate how many of the 25 G = 
l-f2 C IV doublets detected in the Lya forest 
might be false-positive identifications. Two were 
discussed previously. There were 19 from the "def- 
initely C IV" (G = 1) sample. In 15 of the forest 
doublets, the detection of other associated tran- 
sitionj^ and/or multi-component profileJ^ lent 
credibility to the identification as C IV absorption. 

There were eight doublets in the Lya forest 
that had no associated transitions, no Lya cover- 
age, and nondescript profilcs0 Therefore, there 



"The G = 1 21548 = 0.81814 system towards PG1634+706 
is an example of a C IV doublet with associated transitions 
such as Lya and the Si IV doublet. 

^*For example, the G = 1 zi54g = 0.57632 doublet towards 
PGOl 17-1-213 clearly has a multi-component profile, with 
C IV absorption correlated for almost 200kms~^. 

I'^The G = 1 21548 = 0.74843 doublet towards MARK132 has 



36 



Table 10 
SIGHTLINE SELECTION BIAS 



Target 




(kms^i) 


Bias 




(kms^i) 


logAf(C+^) 


PG0117+213 


0.5763" 


600 


Mgll 


1 


3 


14.52 ± 0.03 




0.72907'' 


600 


Mgll 









PG1206+459 


0.93" 


1500 


Mgll 


2 


-503 


> 15.39 












661 


14.25 ± 0.09 


PG1211 + 143 


0.051 


400 


D/H 


1 


39 


14.01 ± 0.03 




0.0652 


400 


D/H 


1 


-228 


13.28 ± 0.05 


PG1241 + 176 


0.5504 


400 


Mgll 


1 


57 


14.60 ± 0.05 




0.5584 


400 


Mgll 


1 


3 


14.26 ± 0.10 




0.8954 


400 


Mgll 


1 


9 


< 13.81 


PG1248+401 


0.7729 


400 


Mgll 


1 


2 


14.66 ± 0.06 




0.8545 


400 


Mgll 


1 


94 


14.61 ± 0.04 


PKS1302-102 


0.0940 


400 


D/H 









CS0873 


0.66 


400 


Mgll 


1 


162 


14.27 ± 0.02 


PG1634+706 


0.701 


400 


D/H, LLS 










IChurchill et al.| (i2000|) 
^Churchill et all lll999bll 

""zabs ^ 0.9254,0.9276, and 0.9342 



Note. — Sightlines targeted for specific absorption systems. We quote the target absorption 
system redshift Zabs from the proposal unless noted otherwise. By default, we excluded ab - 
sorbers with l^i^absl ^ 400 kms""*", the clustering scale measured bv lChurchill fc VogH I I2001D . 
The exceptions are taken from the literature. 



were fewer diagnostics available for evaluating the 
doublet identification. We include these doublets 
in our sample because we estimate the Lya con- 
tamination with Monte-Carlo simulations (§ 14. 3|) . 
There is a total of ten systems in the Lya forest 
that satisfy few of the C IV characteristics (see Ta- 
ble[3]): the eight doublets with nondescript profiles 
and the two that lack significant Lya absorption, 
discussed previously. The simulations of the Lya 
forest contamination rate excludes all ten from be- 
ing Lya forest contamination at the 99.9% confi- 
dence level. 

Ultimately, there were no definite detections of 
C IV doublets without associated Lya absorption. 

7. Summary 

We conducted the largest survey for Z1548 < 1 
C IV absorbers to date. We surveyed 49 sight- 
lines with the HST STIS and/or GHRS archival 
spectra with moderate signal-to-noise ratios and 
resolution. All absorption-line features and actual 
absorption lines were identified by an automated 
feature-finding algorithm. From this list, candi- 



a nondescript profile. 



date C IV systems were assembled, based solely 
on the doublet's characteristic wavelength separa- 
tion (see Table [2]). 

We visually inspected all candidates with rest 
equivalent widths Wr,i548 > 3crvi/^,i548. After con- 
sidering the various diagnostics {e.g., profile; see 
Table [S]), we identified 44 definite (G = 1) and 19 
likely (G = 2) C IV systems. Of these, we only 
analyzed the 38 G = 1 and five G = 2 doublets 
where both lines were detected with Wr > 3crvi/^. 
All subsequent analyses considered the full sample 
and two redshift bins, divided at, approximately, 
the median redshift: Z1548 < 0.6 (20 G = 1-1-2 
doublets) and 0.6 < Z1548 < 1 (23). 

From synthetic spectra, we estimated the un- 
blocked, co-moving pathlength AX to which our 
survey was 95% complete as a function of A^(C"'"^) 
and Wr, 154S (see Figure [7]). For the strongest 
absorbers (logA^(C+3) > 14.1 and Wr,i548 > 
157mA), AA = 15.1. 

There were ten G = 1+2 doublets detected in 
the Lya forest that had nondescript profiles, no 
associated metal lines, and either no Lya absorp- 
tion or Lya coverage did not exist. We estimated 
the contamination rate that Lya forest lines would 
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be mistaken for C IV doublets, with Monte-Carlo 
simulations. The rate is small: dN/dX < 0.11. 
Therefore, it was ruled out at the 99.9% confidence 
level that all of the ten doublets were Lya-as-C IV 
pairs. However, the two G = 2 doublets without 
Lya absorption {Wr^hya < 3(Th/^) were consistent 
with being Lya forest lines masquerading as C IV 
doublets. 

With a maximum-likelihood analysis, we mod- 
eled the column density frequency distribution 
with a power law: f{N{C+^)) = fci4(iVc+3/^o)"" 



= -1.75 



-HO. 28 



0.33 

the 0.6 < 



The best-fit exponent was 

for the ^1548 < 0.6 sample. For _ 
Zi548 < 1 sample, the best- fit value was = 
— I.39I027. There has been no consensus on 
the power-la w exponent in the various zi .^4s > 



1.5 surveys (lEUison et al.' '2000; 'Songailal 12001 



Boksenberg et al . 2003 : .Songaila 2005). 



We measured the C IV absorber line density to 
be dA/biv/dX = 4.lt°] for logiV(C+3) > 13.2 
and Wr,i548 > 52 mA from the full G = H-2 sam- 
ple. dJ\fciv/dX has not evolved significantly since 
z = 5 (see Figure fT4|). 

The sightlines analyzed in the current study 
were observed for a variety of reasons: e.g., the 
targets were UV bright or a known damped Lya 
system (DLA) lay along the sightline. The latter 
reason, and similar ones, might have introduced a 
bias into our survey, since sightlines with DLAs, 
Lyman-limit systems, and Mg II absorbers typi- 
cally correlate wi t h C IV doublets. Assuming that 
iBarlow fc Tvtleil (|l998f) provide an unbiased sur- 
vey for C IV absorbers, we first compared our mea- 
sured redshift densities for 0.2 < Z1548 < 0.8 and 
Wr > 400 mA. We agreed within la and actually 
measured a lower incidence of strong absorbers. 

We also excluded the 11 doublets close to 
the redshift of the targeted systems (see Table 
nop , then re-measured the absorber line density 
and re-fit the frequency distributions. The "un- 
biase d" absorber line density agreed less well 
with iBarlow fc Tvtleil (|l998l l.7a). The best-fit 
f{N{G~^^)) exponent and coefficient differed by 
< la, respectively, compared to the values for the 
full sample. Categorically excluding the C IV dou- 
blets possibly associated with e.g., Mg II absorbers 
biased our survey against strong logiV(C+^) > 14 
doublets. Therefore, we concluded that the results 
for the full C IV sample represented the results of 



a trul y unbiased survey, akin to IBarlow fc: Tvtler 
(|1998[ ). 

Measuring the C+^ mass density relative to the 
critical density at 21548 < 1 was a principle aim 
of this study. High-redshift studies have agreed 
that Qq+s changed little or not at all from z = 
5-^1, once cosmology was properly taken into 
account (Songaila '200lL|Boksenberg et al.' '200: ' 
Schaye et al.. 2003: Pettini et a l. 2003 : Songail; 
2005t IScannapieco et a l. 20061. For the full G = 
1 sample with (z) = 0.654 {tage — 7.5 Gyr), the 
integrated = (6.20^1^2) x 10"* for 13 < 

logA^(C+3) < 15. This was a 2.8 ± 0.7 increase 
over the 1 < z < 5 values. 

We assumed a simple linear model for the tem- 
poral evolution of Qq+s in order to estimate the 
rate of change in the mass density. The linear 
regression for the z < 5 data indicated a slowly 
increasing f2c+3 with increasing age of the Uni- 
verse, at the > 3-a level: dfJc+s /dtage = (0.42 -|- 
0.2) X 10~^Gyr"^ This result relied on the mea- 
surements at z < 1; without which, the slope 
was unconstrained (df^c+a /dtage = (0.15 + 0.3) x 
10-8 Gyr-i). 

Intuitively, it might not appear surprising that 
ric+a continually increases from z = 5 0, since 
the overall metallicity of the Universe increases. 
However, increasing metallicity would not trans- 
late directly to increasing C"'"'^ mass density. The 
C IV mass density is also subject to the changing 
ionizing background in the Universe, which favors 
the C IV transition less at low redshift (see e.g., 
iM adau k Haardtl l2009[ ). lOppenheimer fc Dave 
(2008|) concluded that metallicity and ionizing 
background balanced one another from z « 5 — > 1, 
resulting in a nearly constant flc+s . 

We emphasized that the simple linear model for 
riQ+a over tage was not physically motivated. The 
observed C+^ mass density only accounted for dou- 
blets with 13 < logiV(C+^) < 15, and absorbers in 
this range likely arise in different physical environ- 
ments at high and low redshift. We discussed the 
changing nature of what C IV absorption traces; 
perhaps it probed the low-density IGM at z > 1 
and galaxy halos at z < 1. We will explore, in a 
future paper, the physical conditions and environ- 
ments of C IV absorbers in cosmological simula- 
tions with various feedback prescriptions. 



38 



We thank P. Jonsson for programming ad- 
vice on the Monte-Carlo simulations and the 
maximum- likelihood analysis. We also thank B. 
Oppcnheimer for useful discussions about the 
evolution of f2c+3 and the effects of cosmology. 
Based on observations made with the NASA- 
CNES-CSA Far Ultraviolet Spectroscopic Ex- 
plorer. FUSE is operated for NASA by the Johns 
Hopkins University under NASA contract NAS5- 
32985. Based on observations made with the 
NASA/ESA Hubble Space Telescope Space Tele- 
scope Imaging Spectrograph and Goddard High- 
Resolution Spectrograph, obtained from the data 
archive at the Space Telescope Institute. STScI 
is operated by the association of Universities for 
Research in Astronomy, Inc. under the NASA 
contract NAS 5-26555. The current study was 
funded by the HST archival grant 10679 and the 
NSF CAREER grant AST 05_48180. 

Facilities: FUSE, HST (STIS), HST (GHRS) 

REFERENCES 

Aguirre, A., Hcrnquist, L., Schayc, J., Katz, N., 
Weinberg, D. H., & Gardner, J. 2001, ApJ, 561, 
521 

Bahcall, J. N., et al. 1993, ApJS, 87, 1 

Barlow, T. A., & Tytler, D. 1998, AJ, 115, 1725 

Becker, G. D., Ranch, M., & Sargent, W. L. W. 
2009, ApJ, 698, 1010 

Boksenberg, A., Sargent, W. L. W., & Ranch, M. 
2003, ArXiv Astrophysics e-prints 

Brandt, J. C, et al. 1994, PASP, 106, 890 

Chen, H.-W., Lanzetta, K. M., & Webb, J. K. 
2001, ApJ, 556, 158 

Churchill, C. W., & Le Brun, V. 1998, ApJ, 499, 
677 

Churchill, C. W., Mellon, R. R., Charlton, J. C, 
Jannuzi, B. T., Kirhakos, S., Steidel, C. C, & 
Schneider, D. P. 1999a, ApJ, 519, L43 

— . 2000, ApJ, 543, 577 

Churchill, C. W., Rigby, J. R., Charlton, J. C, & 
Vogt, S. S. 1999b, ApJS, 120, 51 

Churchill, C. W., & Vogt, S. S. 2001, AJ, 122, 679 



Cookscy, K. L. 2009, PhD thesis, Univ. California 

(Santa Cruz) 

Cooksey, K. L., Prochaska, J. X., Chen, H.-W., 
Mulchaey, J. S., & Weiner, B. J. 2008, ApJ, 
676, 262 

Cowie, L. L., Songaila, A., Kim, T.-S., & Hu, 

E. M. 1995, AJ, 109, 1522 

Cui, J., Bechtold, J., Ge, J., & Meyer, D. M. 2005, 

ApJ, 633, 649 

Danforth, C. W., & Shull, J. M. 2008, ApJ, 679, 
194 

Ellison, S. L., Lewis, G. F., Pettini, M., Chaffee, 

F. H., & Irwin, M. J. 1999, ApJ, 520, 456 

Ellison, S. L., Mallen-Ornelas, G., & Sawicki, M. 
2003, ApJ, 589, 709 

Ellison, S. L., Songaila, A., Schaye, J., & Pettini, 
M. 2000, AJ, 120, 1175 

Fernandez-Soto, A., Lanzetta. K. M., Barcons, X., 
Carswell, R. F., Webb, J. K., & Yahil, A. 1996, 
ApJ, 460, L85-I- 

Prye, B. L., Tripp, T. M., Bowen, D. B., Jenk- 
ins, E. B., & Sembach, K. R. 2003, in Astro- 
physics and Space Science Library, Vol. 281, 
The IGM/Galaxy Connection. The Distribu- 
tion of Baryons at z=0, ed. J. L. Rosenberg & 
M. E. Putman, 231 + 

Janknecht, E., Reimers, D., Lopez, S., & Tytler, 
D. 2006, A&A, 458, 427 

Kim, T.-S., Cristiani, S., & D'Odorico, S. 2001, 
A&A, 373, 757 

Lanzetta, K. M., McMahon, R. G., Wolfe, A. M., 
Turnshck, D. A., Hazard, C, & Lu, L. 1991, 

ApJS, 77, 1 

Madau, P., Ferrara, A., & Rees, M. J. 2001, ApJ, 
555, 92 

Madau, P., & Haardt, F. 2009, ApJ, 693, LlOO 

Milutinovic, N., et al. 2007, MNRAS, 382, 1094 

Misawa, T., Tytler, D., lye, M., Storrie-Lombardi, 
L. J., Suzuki, N., & Wolfe, A. M. 2002, AJ, 123, 
1847 



39 



Mobashcr, B. 2002, HST Data Handbook: In- 
troduction to Reducing HST Data, Volume 
1, Version 4.0 (User's Guide, Hubble Space 
Telescope, January 2002. Edited by Baliram 
Mobasher. Baltimore: Space Telescope Science 
Institute) 

Oppenheimer, B. D., & Dave, R. 2006, MNRAS, 
373, 1265 

— . 2008, MNRAS, 387, 577 

Petitjean, R, & Bergeron, J. 1994, A&A, 283, 759 

Pettini, M., Madau, P., Bolte, M., Prochaska, 
J. X., Ellison, S. L., & Fan, X. 2003, ApJ, 594, 
695 

Prochaska, J. X., Chen, H.-W., Howk, J. C, 
Weiner, B. J., & Mulchaey, J. 2004, ApJ, 617, 
718 

Rauch, M., Sargent, W. L. W., Womble, D. S., & 
Barlow, T. A. 1996, ApJ, 467, L5+ 

Ryan-Weber, E. V., Pettini, M., & Madau, P. 
2006, MNRAS, 371, L78 

Ryan-Weber, E. V., Pettini, M., Madau, P., & 
Zych, B. J. 2009, MNRAS, 395, 1476 

Savage, B. D., & Sembach, K. R. 1991, ApJ, 379, 
245 

Scannapieco, E., Ferrara, A., &: Madau, P. 2002, 
ApJ, 574, 590 

Scannapieco, E., Pichon, C, Aracil, B., Petitjean, 
P., Thacker, R. J., Pogosyan, D., Bergeron, J., 
& Couchman, H. M. P. 2006, MNRAS, 365, 615 

Schaye, J., Aguirre, A., Kim, T.-S., Theuns, T., 
Rauch, M., & Sargent, W. L. W. 2003, ApJ, 
596, 768 

Schaye, J., Carswell, R. F., & Kim, T.-S. 2007, 
MNRAS, 379, 1169 

Schaye, J., et al. 2009, ArXiv e-prints 

Simcoe, R. A. 2006, ApJ, 653, 977 

Simcoe, R. A., Sargent, W. L. W., & Rauch, M. 

2004, ApJ, 606, 92 

Songaila, A. 2001, ApJ, 561, L153 



— . 2005, AJ, 130, 1996 
— . 2006, AJ, 131, 24 

Springel, V., & Hernquist, L. 2003, MNRAS, 339, 
289 

Steidel, C. C. 1990, ApJS, 72, 1 

Storrie-Lombardi, L. J., Irwin, M. J., & McMa- 
hon, R. G. 1996, MNRAS, 282, 1330 

Tytler, D., Fan, X.-M., Buries, S., Cottrell, L., 
Davis, C., Kirkman, D., & Zuo, L. 1995, in QSO 
Absorption Lines, ed. G. Meylan, 289 — h 

Wise, J. H., & Abel, T. 2008, ApJ, 685, 40 



This 2-column preprint was prepared with the AAS lATJ^X 
macros v5.2. 



40 



A. Velocity Plots 
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Fig. 17. — Velocity plots of G = 1 C IV systems. The regions of spectra around each absorption line are 
aligned in velocity space with respect to the rest wavelength of the transition and Z1548 • Saturated transitions 
are indicated with the (red) 'S'; transitions with iV(C+'') < 3crjv(c+3) are indicated with the (red) 'W.' The 
regions used to measure Wr and log N are shown by the dark outline. The flux at zero and unity are shown 
with the dash-dot lines (blue and green, respectively); the (black) vertical dashed line indicates v = Okms~^. 



41 



PKS0405-123 (04:07:48.42) 2,548 = 0.36071 PKS0454-22 (04:56:08.91) 2,548 = 0.24010 
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. 17. — G = 1 velocity plots (continued) 
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. 17. — G = 1 velocity plots (continued) 
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17. — G = 1 velocity plots (continued) 
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Fig. 17. — G = 1 velocity plots (continued) 
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Fig. 17. — G = 1 velocity plots (continued) 
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PG1241+176 (12:44:10.83) LyaForest z,54g = 0.55070 PG1241+176 (12:44:10.83) LyaForest 2,548 = 0.55842 
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. 17. — G = 1 velocity plots (continued) 
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Fig. 17. — G = 1 velocity plots (continued) 
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17. — G = 1 velocity plots (continued) 
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CS0873 (13:19:56.28) z,54s = 0.73385 PG1630+377 (16:32:01.10) LyaForest z.jjj = 0.91440 
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. 17. — G = 1 velocity plots (continued) 
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Fig. 17. — G = 1 velocity plots (continued) 
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H1821+643 (18:21:57.10) 2,548 = 0.22503 PHLlSll (21:55:01.53) 2,548 = 0.08091 
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Fig. 17. — G = 1 velocity plots (continued) 
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PG0117+213 (01:20:17.25) LyaForest z,54g = 0.51964 PKS0232-04 (02:35:07.30) LyaForest 2,548 = 0.86818 
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. 18. — Velocity plots of G = 2 C IV systems. (See Figure flTl for description of velocity plot.) 



53 



HE0515-4414 (05:17:07.53) 

1.2 - 



0.5 - 



■o 



1.2 - 



o 

2; 



0.5 



LyaForest z,54g = 0.50601 HE0515-4414 (05:17:07.53) 

1.2 



Lya Forest 



Z,54s = 0.73082 



CIV 1548' 




0.5 



ft. 



S 1.2 
5 

2; 



0.5 





1 1 

1 


1 


I 

CIV 1548 ; 






^ ! ' 

1 


CIV 1550" 




1 ■ 1 



-100 100 

Relative Velocity (km s^') 

HS0747+4259 (07:50:54.63) LyaForest z,;4j = 0.83662 HS0810+2554 (08:13:31.32) 



100 100 

Relative Velocity (km s'') 



Lya Forest 



z,54j = 0.87687 




-100 100 

Relative Velocity (km s ') 



-100 100 

Relative Velocity (km s'') 



18. — G = 2 velocity plots (continued) 
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MARK132 (10:01:29.76) LyaForest z,54g = 0.70776 3C249.1 (11:04:13.77) LyaForest 2,548 = 0.02616 
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. 18. — G = 2 velocity plots (continued) 
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PG1241+176 (12:44:10.83) LyaForest z,54g = 0.48472 PG1248+401 (12:50:48.32) LyaForest 2,548 = 0.55277 
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Fig. 18. — G = 2 velocity plots (continued) 
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Fig. 18. — G = 2 velocity plots (continued) 
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B. MELximum-Likelihood Analysis 

We used the maximum liiceliiiood metliod to fit tlie frequency distributions with a power-law function. 
For this discussion, we will focus primarily on deriving the maximum-likelihood function for the column 
densities, where one must account for lower limits due to saturation effects. Saturation obviously does not 
occur for the equivalent width measurements, and the maximum-likelihood function L derived here r educe s 
to a formalism applicable to W^. The following derivation draws heavily from lStorrie-Lombardi et al.l ( 19961 ). 

Assuming Poisson counting statistics, the probability of detecting M(Ni) absorbers with column densities 
Ni is determined by the Poisson probability distribution: 

P(A/-(7V.);M») = e"'^'^^, (Bl) 

where [li is the expected number, based on the parent distribution. The likelihood function is defined as the 
product of the probability of observing each absorber with 

i i 

Assume that the expected number [li depends on the frequency distribution: 

/i, = /(7V,)AA(A,)AAf, (B3) 

(see Equation [TT|) . Then, consider the limit where the "volume" AA(Ai) AA^ contains at most one absorber. 
With M ^ p (single-detection) +g (no-detection) volumes: 



ln£ = -/ /(A,)AA(AOdA, + Vln(/(A,)AA(A,)), (B4) 



in the limit where AA^ (AA(Ai) is essentially a weight and cannot shrink to zero for all Ni). 

However, our absorbers include column density measurements where we have only a lower limit on 
Ni > Nsat- We need a probability for detecting q saturated absorbers above the saturation limit 
^sat = 10^'^ '^ cm~^, P{q;fi^). Thus, the likelihood function becomes: 

The mean number of saturated absorbers is expected to be: 

A4at - / f{N,)AX{N,)dN, (B6) 

Let = A^sat, the expected number of saturated absorbers from the integral of f{Ni) between Nsat and 
infinity. Thus, taking the natural logarithm and expanding Equation IB5I reduces to: 

\nC = - I f{N,)AX{N,)dN,~ I f{N)AX{N,)dN, (B7) 

{N)dN 



N^i„ JN 



P / /•oc 

Vln(/(A,)AA(7V,))+9ln / /(^^^^ 



Substituting the power-law form of /(A^(C+'^)) (see Equation [T^ provides the maximum-likelihood function 
used in the current study. 
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C. Adjusting Mass Density 

In order to compare 0^+3 from the current study to other studies, we had to account for differences in 
cosmology and range of N{C^^) included in the measurement. Changing the adopted Hubble constant Hq is 
a simple matter of scaling 51^+3 by the ratio of the old Hq to the new (see Equation[Tni recall that pc,o oc Hq 
and AX is independent of Hq). Less simple is adjusting for changes in and JIa, since that enters the 
computation of ^Iq+s in the estimate of the co-moving pathlength AX (see Equation [9|) . Previous authors 
have frequently scaled the summed flc+3 by the ratio of the new pathlength to the old as follows: 



O^j--, -ffo.old— ;~((^)' ^M.old, ^^A,old) 



(CI) 



-ffo,now~; — ^M,ncwi ^^A,ncw) 

dz 



where dX/dz is the derivative of Equation [5] evaluated at the median redshift (z) and with the appropriate 
cosmology. 

This adjustment is an approximation because AX is the sum of parts of spectra that satisfy the redshift 
and column density constraints (see Equation II 0[). Therefore, to correctly adjust for cosmology, each snippet 
SX must be calculated in the desired cosmology and then added together to re-estimate AX. Most published 
studies do not provide the necessary information to do this, which is why the above equation is the standard 
adjustment. However, the approximation works well when the redshift bin is small. For values quoted in the 
text and Figures fT3l and \TE\ we use Equation IC II to adjust for differences in cosmology. 

For m ost of the studie s, th e adjustment i s small since their adopted cosmology and ours were similar. The 



results of lSongailal (j200ll ) and lPettini et al.l (|2003i) were decreased by w 55%, since they adopted an Einstein- 



de Sitter cosmology (EdS) with Hq = 65 kms"""^ Mpc~^, Qm = 1, and ~ 0. Since we cannot correctly 
compute their results precisely for our adopted cosmology, we calculated our rj^+s in their cosmology, by re- 
constructing AX(A^(C+'^)) and re- calculating the summed fic+a. The difference between the values quoted 
in Table [6] and our EdS adjusted by Equation ICll is < 8% for the full sample. Thus, Equation ICll 

adequately adjusts for differences in cosmology when applied to summed values of f2c+3 

The column density limits used in the summed f2c+3 estimation was the other factor we considered 
when adjusting the values from other studies. As mentioned in § 15.31 the high column density doublets 
contain the majority of the mass and dominate the C"''^ mass density. Since there has been no break in the 
frequency distribu t fon, w e limited the integrated measurement to 13 < logA(C"'"'^) < 15. However, 

Scannapieco et all (|2006l ) estimated f7c+3 = (7.54 ± 2.16) x 10"^ by summing their observed doublets with 



12 < log7V(C"'"^) < 16. To fairly compare their result with the current survey, their measured value was 
decreased by 45%. This factor was determined by assuming their summed £7^+3 scaled as the integrated 
f2c+3 : 

''C+3,now -"max, now -"mm, new 



f^c+3,oid Ait"^„ - TV^t"! 



(C2) 



max, old mm, old 

where the new column density limits are 13 < log A(C+'^) < 1 5, to match the cur rent work, and aA? is the 
value measured by the "old" study (also see Equation Wh. F or Scannapieco et al. 120061, the best-fi t gjy 



1. 8. We also made ad j ustm ents (< 30%) to lSongailal (j200l[) . |Pettini et all (|20mlBoksenberg et al.l (|2003l ). 



and 'Pa nforth fc Shulll (2008") since iVmm < 10"^^ cm ^ in their studies and, for most, A^max < 10^^ cm ^, 
which increases the significance of the lower column density absorbers. 

We did not extrapolate these studies and change Nmax ■ These surveys were surely sensitive to doublets 
with log A(C+'^) > 14, but they did not survey a sufficiently large pathlength to encounter them, since they 
are rare. The affect of pathlength is accounted for in the summed il^+s f Equation 1^0)) . where it essentially 
weights the sum of the column densities. Therefore, for these studies, the summed flc+3 values occasionally 
reflect the fact that the strongest absorbers are rare, though they dominate the C"'"'^ mass density. 



'Using Equation [CT] to adjust our EdS, integrated f!Q+3 values introduces errors up to 25%. 
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The result from Ryan- Weber et all ( 20091 ) could have been increased by a factor of 1.4, assuming ajv = 
—1.8 (less for a shallower power law); they were only sensitive to doublets with 13.8 < log7V(C"'"'^) < 15. 
However, we did not adjust their value since we did not know the appropriate slope and did not include the 
highest redshift measurements in the linear regressions to Qq+s over tage- 

The total adjustment for cosmo logy and/or c olumn density hmits for the z > 1 studies cited in thi s paper 
Songailal feoOl): 38%. IPettini et all (l2003l): 81%. iBoksenberg et al.l (|2003[) : 45%, 



are as follows: 



40%, 



Scannapieco et all 92%. lOanforth fc Shulll (|2008h : and 105%. lBecker et al.l 1^200^ . 



D. Detailed Comparison with iMilutinovic et aU (|2007h and iDanforth Sz Shulll (|2008h 

Milutinovic et all (|2007l ) and lPanforth &: Shulll l|2008h surveyed archival STIS E230M and E140M spectra, 



respectively, for intergalactic absorption lines. It is useful, therefore, to compare their results with our own 
as a consistency check. 



In general, our search algorithms and sample selection were consistent with IMilutinovic et al.l (|2007f) and 
IDanforth &: Shull (2008), who used different procedures. The current survey, however, has the advantage of 
analyzing the largest sample of Z1548 < 1 sightlines in a consistent fashion. 



D.l. IMilutinovic et al.l (|2007h 



Milutinovic et al.l (|2007t ) focused on eight sightlines with STIS E230M observations and compiled compre- 
hensive line lists for each sightline. There were 24 C IV systems in their survey. We agreed with 19 of these, 
including the associated transitions {e.g., Lya, Si II 1260), though we did not search for all of the transitions 
that they did {e.g., C II 1334, Fe I I 1608). We coun t ed th e DLA Z1548 = 0.92677 C IV doublet towards 
PG1206+459 as one system, while Milutinovic et all (2007) divide d it into two. Th r ee of the remaining 
five systems are acknowledged to be questionable identifications bv IMilutinovic et "aL ( 20071 ) and were not 
included in our G — 1+2 group at all. Of the other two, the reported doublet at 21543 = 0.7760 towards 
PG1248+401 does not show up in our survey because the 1548 line was blended and the 1550 would not be 
detected with Wr, 1550 > ScrWr-- The other unconfirmed doublet at 21548 = 0.9903 towards PG1634+706 was 
associated with the Mg II absorber, which was the motive for observing the sightline (see i) l6.4[) . This doublet 
would lie in the highest-wavelength order of the E230M spectrum, which we excluded due to questionable 
data quality flags. 

Considering the opp osite comparison, we identified 29 C IV doublets in the eight sightline s that 
Milutinovic et al. ( 20071 ) also analyzed. Of these, 11 were identified as other transitions bv [Milutinovic et al.l 
(|2007l ). In nine cases, the absorption lines we list as C IV doublets were either: listed as tentative Lya 
lines; not identified; or not listed at all. The two remaining disputed cases are the ones we identify as 
Z1548 = 0.48 4 72 towards P G1241-H76 and z, m^ = 0.55277 towards PG1248-h401, both G = 2 doublets. 



For the PG1241-M76 lines, IMilutinovic et al.1 (|2007l ) identifies our 1548 hue as O VI 1037 a,t z = 1.2142 
and 1550 as a tentative Lya line. The 1548 profile is c onsistent with being blen ded, so we conclude that 
both identifications are valid. For the PG1248-I-401 lines. IMilutinovic et al.l (|2007l ) lists our 1548 line as Lya 
absorption, with no associated transitions, and does not list our 1550 line, which we detect at > SaWr-- We 
stand by our identification and accept the possibility that the 21543 = 0.55277 doublet is Lya contamination, 
as discussed in S 14.31 



D.2. IDanforth &: Shulll (|2008h 

Danforth fc Shulll (|2008f ) focused on 28 sighthnes with STIS E140M observations and supplementary FUSE 



et al. 1 1 I2OO3D was a conference proceeding, which we do not discuss in detail. To summarize, they surveyed nine sightlines 
with STIS E140M spectra. We agreed with seven of their nine identified C IV doublets. Though we agreed well with their 
dA/c IV /dz over equivalent width limit, their reported 0(-.+3 value was almost double our value for the 2 < 0.6 sample. However, 
this is likely due to small number statistics, since we agreed (< Icr) when we focused just on the E140M data (see Table|6ll. 
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spectra. They conducted a Lya-targeted survey, where they identified Lya absorbers first, then searched 
for other, associated intergalactic transitions {e.g., Lyf3, C III). They detected 24 Z1548 < 0.12 C IV doublets 
with at least one line detected at > Aaw,. so long as Rw was consistent (see § 13. 2p . For a total pathlength 
Az — 2.42, they measured dA/civ/dz — lO^j ^'^^ > 30 mA. They also fit a power law to dA/civ/dz 
in bins of AlogiV(C+3) = 0.2; the best-fit exponent was a AT = —1.79 ±0.17. They measured a summed 
= (7.78 ± 1.47) X 10-** for 12.83 < log7V(C+3) < 14.13. 

The current survey included the same 28 sightlines th ough not all of the FUSE data, since some had 
S/N too low for continuum fitting. Of the 24 doublets in IPanforth fc ShuU (2008), we detected nine also 
and could neither agree or disagree with three, which were detected in the zeroth order of the STIS E140M 
spectra, which we did not include due to questionable data quality flags. The zeroth order covered 1711 A < 
A < 1729 A (or 0.105 < z < 0.115), and the sensitivity of the E140M spectra drops rapidly at A > 1700 A. 
Therefore, we would not likely have inc luded any doublet s in th e zeroth order in our analysis, since both 
lines were not detected with Wr > 'Sa in iDanforth &: Shulll (|2008h . 

Of the remaining 11 doublets, two towards NGC7469 were excluded from our sample because they were 
within 3000 km s~^ of the background source. The other nine either did not look like absorption lines or 
th e C IV 1548 line did n ot have Wr.iR48 > 3gvi/ in our co-added and normalized spectra. As discussed 



Cooksev et all (|2008h . IPanforth fc Shulll (j2008l ) cited questionably small errors (< 5%) for their rest 



equivalent widths, which led them to consider more lines to be detected at a given significance level. 

We detected 12 C IV doublets in the E140M spectra, that includes three doublets not cited in 
IDanforth fc Shulll (l2008l) . For two cases, we identified the absorption lines as G = 2 C IV doublets with 
only the 1548 line detected with Wr,i548 > 3(Hv,: Z1548 = 0.02616 towards 3C249.]E1 and Z1548 = 0.08077 
towards HS1700-I-6416. These doublets were not included in our analyses. The last case was the = 1, 
•^1548 = 0.00574 doublet towards QSO-123050-1-011522, which was also associated with a Si IV doublet (an- 
other target absorption lines of their study). They did not detect this doublet because the Ly/3 line was 
contaminated by a spurious artifact in the FUSE spectra, but they now agree with our identification (C. 
Danforth, private communication) 



^ ^Danforth fc Shulll ||2008| ') identified the > Saw,, absorption line (that wc list as C IV 1548) as Lya at z = 0.30788, which was 
their default identification for single lines. 
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